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Abstract 

In this paper, we show how to reduce the Penrose conjecture to the 
known Riemannian Penrose inequality case whenever certain geometri- 
cally motivated systems of equations can be solved. Whether or not these 
special systems of equations have general existence theories is therefore 
an important open problem. The key tool in our method is the deriva- 
tion of a new identity which we call the generalized Schoen-Yau identity, 
which is of independent interest. Using a generalized Jang equation, we 
propose canonical embeddings of Cauchy data into corresponding static 
spacetimes. In addition, our techniques suggest a more general Penrose 
conjecture and generalized notions of apparent horizons and trapped sur- 
faces, which are also of independent interest. 



1 Introduction 

In addition to their intrinsic geometric appeal, the Penrose conjecture [53] and 
the positive mass theorem [21] are fundamental tests of general relativity as a 
physical theory. In physical terms, the positive mass theorem states that the 
total mass of a spacetime with nonnegative energy density is also nonnegative. 
The Penrose conjecture, on the other hand, conjectures that the total mass of 
a spacetime with nonnegative energy density is at least the mass contributed 
by the black holes in the spacetime. In this section, we will explain how these 
simple physical motivations translate into beautiful geometric statements. 

After special relativity, Einstein sought to explain gravity as a consequence of 
the curvature of spacetime caused by matter. In contrast to Newtonian physics, 
gravity is not a force but instead is simply an effect of this curvature. As an 
analogy, consider a heavy bowling ball placed on a bed which causes a significant 
dimple in the bed. Now roll a small golf ball off to one side of the bowling ball. 
Note that the path of the golf ball curves around the bowling ball because of the 
curvature of the surface of the bed. In this analogy, the bowling ball represents 
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the sun, the golf ball represents the earth, and the surface of the bed represents 
spacetime. Whereas Newton explained the curvature of the path of the smaller 
object by asserting an inverse square law force of attraction between the two 
objects, Einstein declared that the curvature of the smaller object's trajectory 
was due to the curvature of spacetime itself, and that objects which did not 
have forces (other than gravity) acting upon them followed geodesies in the 
spacetime. That is, according to general relativity, the sun and all of the planets 
are actually following geodesies, curves with zero curvature, in the spacetime. 

It should also be noted that general relativity is entirely consistent with 
large scale experiments, whereas Newtonian physics is not. The most notable 
example may be the precession of the orbit of Mercury around the Sun. Whereas 
general relativity predicts the rate at which the elliptical orbit precesses around 
the Sun to as many digits as can be measured, Newtonian physics is off by 
almost one percent, with all possible excuses for the discrepancy having been 
eliminated. The question, then, is how to turn the beautiful and experimentally 
verified idea of matter causing curvature of spacetime, which Einstein called his 
happiest thought, into a precise mathematical theory. 

First, assume that {N*,gN) is a Lorentzian manifold, meaning that the 

metric has signature ( h ++) at each point. Note that at each point, 

time-like vectors (vectors v with gpf{v,v) < 0) are split into two connected 
components, one of which we will call future directed time-like vectors, and the 
other of which we will call past directed time-like vectors. 

Next, define T{v, w) to be the energy density going in the direction of v 
as measured by an observer going in the direction of w, where v, w are future- 
directed unit time-like vectors at some point p £ N. In addition, suppose that 
T is linear in both slots so that T is a tensor. Then the physical statement that 
all observed energy densities are nonnegative translates into 

T{v,w) > 

for all future-directed (or both past-directed) time-like vectors v and w at all 
points p £ N, known as the dominant energy condition. 

The goal, then, is to set T, which is called the stress-energy tensor, equal to 
some curvature tensor. A natural first idea is to consider the Ricci curvature 
tensor since it is also a covariant 2-tcnsor. In fact, this was Einstein's first idea. 
However, the second Bianchi identity on a manifold N with metric tensor g^ 
implies that 

div(G) = 0, 

where G = Ricn — ^RnSn, RicN is the Ricci curvature tensor, and Rn is the 
scalar curvature. This geometric identity led Einstein to propose 

G = SttT, (1) 

known as the Einstein equation, since as an added bonus we automatically get 
a conservation-type property for T, namely div(r) = 0. Naturally this is a very 
nice feature of the theory since energy and momentum (the spatial components 
of the energy vector) are conserved in every day experience. 
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The next step in pursuing this hne of thought is to try to find examples 
of spacetimes which satisfy the dominant energy condition, the simplest case 
of which would be spacetimes with G = which are naturally called vacuum 
spacetimes. Taking the trace implies that such spacetimes (in 2+1 dimensions 
and higher) have zero scalar curvature and therefore zero Ricci curvature as 
well. The first example (in 3+1 dimensions) is clearly Minkowski space 

(R^, -~dt^ + dx^ + dy^ + dz^) 

which has zero Riemann curvature tensor. The second simplest example of a 
spacetime with G = 0, 

j^R X (R-^ \ i?,„/2(0)), - ' '^^^ + + ^'^''^ + '^^^ + '^'^^ ' 

(2) 

where r = 1/2;^+?/^ + 2^, is a one parameter family of spacetimes called the 
Schwarzschild spacetimes. When m > 0, these spacetimes represent static black 
holes in a vacuum spacetime. 

While the Schwarzschild spacetime can be covered by a single coordinate 
chart (see Kruskal coordinates described in section [21) , the coordinate chart 
above only covers the exterior region of the black hole and has a coordinate sin- 
gularity (not an actual metric singularity) on the coordinate cylinder r = m/2. 
For our purposes, however, wc will only be interested in the exterior region of 
the Schwarzschild spacetime, which physically corresponds to the region where 
observers have yet to pass into the event horizon of the black hole, which is the 
point of no return from which not even light can escape back out to infinity. 

Spacetimes which may be expressed in the form 

(R X M, ~(\)(xfdt^ + .g) , 

where t e R, .t G M , and g is a Riemannian (positive definite) metric on M, 
are called static spacetimes. This name is appropriate since we see that the 
components of the spacetime metric in this coordinate chart do not depend on t 
but instead are entirely functions of x. Note also that static metrics are defined 
not to have any time/spatial cross terms. (Spacetimes which allow time/spatial 
cross terms but where the metric components still only depend on x are called 
stationary spacetimes.) 

An important result, first proved by Bunting and Masood-ul-Alam [6] us- 
ing a very clever argument involving the positive mass theorem, is that the 
only complete, asymptotically flat static vacuum spacetimes with black hole 
boundaries (or no boundary) are the two spacetimes that we have listed so far, 
Minkowski and Schwarzschild. This fact suggests that a thorough understand- 
ing of these two spacetimes, including what makes them special as compared 
to generic spacetimes, may be important for understanding some of the most 
fundamental properties of general relativity. 

In fact, the Minkowski and Schwarzschild spacetimes are the extremal space- 
times for the positive mass theorem and the Penrose conjecture, respectively. 
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That is, the case of equahty of the positive mass theorem states that any space- 
hke hypersurface of a spacetime satisfying the hypotheses of the positive mass 
theorem which has m = can be isometricahy embedded into the Minkowski 
spacetime. Similarly, the case of equality of the Penrose conjecture (which, while 
still a conjecture, has no known counter-examples in spite of much examination) 
states that any space-like hypersurface of a spacetime satisfying the hypotheses 
of the Penrose conjecture which has m = ^ AjX^-n (or to be more precise, the 
region outside of the outermost minimal area enclosure of the apparent horizons) 
can be isometrically embedded into the Schwarzschild spacetime. 

Before we can state these theorems, though, we need to define a few terms. 
The basic object of interest in this paper is a space-like hypersurface of a 
spacetime N"^^ along with the induced metric g on and its second funda- 
mental form fc in the spacetime. 

From this point on we will assume that Ivp has a global future directed 
unit normal vector n future in the spacetime. This standard assumption is not 
necessary for stating the dominant energy condition or, as we will see in section 
[31 for defining generalized apparent horizons or stating the generalized Penrose 
conjecture, but it is necessary for the traditional definition for apparent horizons 
of black holes, as wc will sec. So for convenience, wc will abuse terminology 
slightly and also call 

k{V,W) = -{VvW, 

'^future) (3) 

the second fundamental form of , where V, W arc any vector fields tangent 
to and V is the Levi-Civita connection on the spacetime N"^. In this manner 
we are defining fc to be a real-valued symmetric 2-tensor, where the true second 
fundamental form, which takes values in the normal bundle to M^, is k-n future ■ 

Definition 1 The triple {M^,g,k) is called the Cauchy data of for any 
positive definite metric g and any symmetric 2-tensor k. 

This name is appropriate because this is the data required to pose initial 
value problems for p.d.e.'s such as the vacuum Einstein equation G = 0, or the 
Einstein equation coupled with equations which describe how the matter evolves 
in the spacetime. Note when is flowed at unit speed orthogonally into the 
future that 

|,,«^».,„ 

so that k is in fact the first derivative of g in the time direction (up to a factor). 

Curiously, as we will see, the positive mass theorem and the Penrose con- 
jecture reduce to and are fundamentally statements about the Cauchy data of 
space-like hypersurfaces of spacetimes, not the spacetimes themselves. 

Definition 2 At each point on M^, define /i = Tin futures ^future) to he the 
energy density and the covector J on to be the momentum density, where 
J{v) = Tlufuture^v), where v is any vector tangent to M. 
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By the Einstein equation (equation[T|) and the Gauss-Codazzi identities [55], 
it follows that /x and J can be computed entirely in terms of the Cauchy data 
{M^,g,k). In fact, 

(Stt) /I G(n/„t„^e,n/uture) = (i? + tr (fc)^ - 1 1 fc f )/2 (4) 
(Stt) J = G{nfuture,-) = div (fc - tr(fc)g) , (5) 

where R is the scalar curvature of {M^,g) at each point, and the above traces, 
norms, and divergences are naturally taken with respect to g and the Levi-Civita 
connection of g. Then the dominant energy condition on T implies that wc must 
have 

f^>\Jl (6) 

which wc will call the nonncgative energy density condition on (M^ , 9,^), where 
again the norm is taken with respect to the metric g on M'^ . 

Equations [3] and [5] are called the constraint equations because they impose 
constraints on the Cauchy data {M^,g,k) for each initial value problem. For 
example, we clearly need to impose /i = and J = on any Cauchy data which 
is meant to serve as initial conditions for solving the vacuum Einstein equation 
G ~ 0. However, for our purposes throughout the rest of this paper, we will be 
interested in Cauchy data {M^, g, k) which only needs to satisfy the nonncgative 
energy density condition in inequality [51 Since the assumption of nonncgative 
energy density everywhere is a very common assumption, the theorems we prove 
will apply in a very broad set of circumstances. 

Next we turn our attention to the definition of the total mass of a spacctimc. 
Looking back at the Schwarzschild spacetimc, time-like geodesies (which rep- 
resent test particles) curve in the coordinate chart as if they were accelerating 
towards the center of the spacetime at a rate asymptotic to vnjr^ in the limit 
as r goes to infinity. Hence, to be compatible with Newtonian physics (with the 
universal gravitational constant set to 1) in the low field limit, we must define 
m to be the total mass of the Schwarzschild spacetime. 

More generally, consider any spacetime which is isometric to the Schwarzschild 
spacetime with total mass m for r > rg and which is any smooth Lorcntzian 
metric satisfying the dominant energy condition on the interior region. Of course 
the Schwarzschild spacetime satisfies the dominant energy condition since it has 
(7 = 0. Then the same argument as in the previous paragraph applies to this 
spacetime, so its total mass must be m as well. This last example inspires the fol- 
lowing definition, which comes from considering the t = slice of Schwarzschild 
spacetimes. 

Definition 3 The Cauchy data {M'^,g,k) will be said to be Schwarzschild at 
infinity if can be written as the disjoint union of a compact set K and a 
finite number of regions Ei ( called ends ), where k — on each end and each 

{Ei,g) is isometric to (rJ^ \ Br. (0), (l 4- [dx'^ + dy^ + dz'^fj for some rui 

and some Ri > max(0, —771^/2). In addition, the mass of the end Ei will be 
defined to be rrii. 
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We refer the reader to [17] and [52] for more general definitions of asymp- 
totically flat Cauchy data, but for this paper the special case of being precisely 
Schwarzschild at infinity is sufficiently interesting. 

Typically we will be interested in Cauchy data with only one end. However, 
sometimes it is convenient to allow for the possibility of multiple ends. Each 
end represents what we would normally think of as a spatial slice of a universe, 
and the positive mass theorem and the Penrose conjecture may be applied to 
each end independently. In fact, since ends can be compactified by adding a 
point at infinity and then using a very large spherical metric on the end without 
violating the nonnegative energy density assumption, without loss of generality 
we may assume that any given Cauchy data has only one end for the problems 
we will be considering. 

Theorem 1 (The Positive Mass Theorem, Schoen-Yau, 1981 fF^ ; Witten, 
1981 JMI) 

Suppose that the Cauchy data {M^^g,k) is complete, satisfies the nonnegative 
energy density condition n > \J\, and is Schwarzschild at infinity with total 
mass m. Then 

m > 0, 

and m = if and only if {M^,g, k) is the pullback of the Cauchy data induced 
on the image of a space-like embedding of into the Minkowski spacetime. 

The above theorem has an important special case when fc = which is al- 
ready extremely interesting. Note that the nonnegative energy condition reduces 
to simply requiring (Ajf^,^) to have nonnegative scalar curvature. 

Theorem 2 (The Riemannian Positive Mass Theorem, Schoen-Yau, 1979 fii6] /; 
Witten, 1981 JS^) 

Suppose that the Riemannian manifold {M^,g) is complete, has nonnegative 
scalar curvature, and is Schwarzschild at infinity with total mass m. Then 

m > 0, 

and TO = if and only if {M^, g) is isometric to the flat metric on R'^. 

The adjective Riemannian was introduced by Huiskcn-Ilmanen in [14j since 
the theorem is a statement about Riemannian manifolds as opposed to Cauchy 
data in the more general case. We remind the reader that Cauchy data (M^, g, k) 
is still required to have a positive definite metric g. 

Notice that the Riemannian positive mass theorem is a beautiful geometric 
statement about manifolds with nonnegative scalar curvature. In fact, Schoen- 
Yau were studying such manifolds [29| for purely geometric reasons when they 
first realized that they could use minimal surface techniques to prove the Rie- 
mannian positive mass theorem. They then observed [17] that theorem[T] (which 
is quite mysterious from a geometric point of view without the physical moti- 
vation) reduced to theorem [2] after solving a certain elliptic p.d.e. on [M^ ,g, h) 
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called the Jang equation, named after the theoretical physicist who first intro- 
duced the equation in |15j . 

Witten's proof of the positive mass theorem uses spinors and proves both 
of the above statements by applying the Lichnerowicz-Weitzenbock formula to 
a spinor which solves the Dirac equation, and then integrating by parts. This 
proof has a strong appeal because it computes the total mass as an integral of a 
nonnegative integrand. However, so far it has not been clear how to generalize 
this approach to achieve the Penrose conjecture, although very interesting works 
in this direction include [TT] and [TO] . 

Before we can state the Penrose conjecture, we need several more definitions. 
For convenience, we modify the topology of by compactifying all of the ends 
of except for one chosen end by adding the points {ook}- (However, the 
metric will still not be defined on these new points.) 

Definition 4 Define S to be the collection of surfaces which are smooth compact 
boundaries of open sets U in , where U contains the points {oo^.} and is 
bounded in the chosen end. 

All of the surfaces that we will be dealing with in this paper will be in 
S. Also, we see that all of the surfaces in S divide into two regions, an 
inside (the open set) and an outside (the complement of the open set). Thus, 
the notion of one surface in S enclosing another surface in S is well defined as 
meaning that the one open set contains the other. 

Definition 5 Given any S G 5, define T, £ S to be the outermost minimal area 
enclosure of S. 

That is, in the case that there is more than one minimal area enclosure of the 
surface S, choose the outermost one which encloses all of the others. The fact 
that an outermost minimal area enclosure exists and is unique roughly follows 
from the following: if dA and dB are both minimal area enclosures of some 
surface, then so are d{A U B) and d{A n B) since \d{A D B)\ + \d{A n B)\ = 
\dA\ + \dB\ ~ 2Amin and both have area at least Amm- A rigorous proof that 
the outermost minimal area enclosure of a surface in an asymptotically flat 
manifold exists and is unique is given in jl4j . 

Definition 6 Define S G 5 m (M^,g,k) to be an apparent horizon if it is one 
of the following three types of horizons, 

a future apparent horizon if 



i^s+trs(fc) =0 



on E, 



(7) 



a past apparent horizon if 



-ffE -trs(fc) =0 



071 S, 



(8) 



and a future and past apparent horizon if 



iJs = and tr5](fc) =0 onT,, 



(9) 
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where iJs *s the mean curvature of the surface E in (AI^ , g) ( with the sign 
chosen to be positive for a round sphere in flat 'Rp) and trs{k) is the trace of k 
restricted to the surface E. 

Note that equation [9] follows from assuming both equations [7] and [8] everywhere 
on E. Also note that E is not required to be connected, although from a physical 
point of view each component of E is usually thought of as the apparent horizon 
of a separate black hole. Finally, observe that all three types of horizons are 
simply minimal surfaces (surfaces with zero mean curvature) in the important 
special case when fc = 0. 

Physically, the only relevant apparent horizons for a spacecraft flying around 
in a spacetime are future apparent horizons, because spacecraft are only con- 
cerned about being trapped inside black holes in the future. Mathematically, 
however, merely changing the choice of global normal vector n future to in 
N'^ to ~n future changes the sign on k which causes past apparent horizons to 
become future apparent horizons, and vice versa. 

Equations El [H [9] are actually all conditions on the mean curvature vector 
of E in the spacetime. Note that at each point of E^, the normal bundle, of 
which the mean curvature vector is a section, is a 2-dimensional vector space 
with signature ( — h). Naturally, a basis for this vector space is any outward 
future null vector along with any outward past null vector. Since Y? bounds 
a region in , outward is well-defined, and since there exists a global normal 
vector n future to M^, the future direction is well-defined. 

Geometrically, if one fiows a submanifold in the normal directions fj, then 
the rate of change of the area form of the submanifold is given by 

j^dA=-{rl,H)dA 

where H is the mean curvature vector. It turns out that the mean curvature of 
a surface E^ contained in a slice with Cauchy data {M^,g,k) has coordinates 
(trs(fc), — iJ), where the first coordinate is in the unit future normal direction 
to the slice and the second component is in the unit direction outward perpen- 
dicular to the surface and tangent to the slice. With this convention, then the 
vector with components (1,1) is an outward future null vector, and the vector 
with components (-1,1) is an outward past null vector. 

Hence, equation [7] is equivalent to requiring that, at each point of E, the 
dot product of the mean curvature vector with any outward future null vector 
is zero (which implies that the mean curvature vector is a real multiple of the 
outward future null direction). Similarly, equation [8] is equivalent to saying that 
the dot product of the mean curvature vector with any outward past null vector 
is zero. Thus, future apparent horizons have the property that their areas do not 
change to first order when flowed in outward future null directions. The same 
is true for past apparent horizons when flowed in outward past null directions. 

We are now able to state the Penrose conjecture. An excellent survey of this 
conjecture is found in [^D] . 
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Conjecture 1 (The Penrose Conjecture, 1973 \25^ - Standard Version) 
Suppose that the Cauchy data {M^,g,k) is complete, satisfies the nonnegative 
energy density condition fi > \J\, and is Schwarzschild at infinity with total 
mass m in a chosen end. IfYp'ES is a future apparent horizon, then 

where A is the area of the outermost minimal area enclosure = dU^ of . 
Furthermore, equality occurs if and only if {M^ \ U^, g, k) is the pullback of the 
Cauchy data induced on the image of a space-like embedding of \ into 
the exterior region of a Schwarzschild spacetime (which maps to a future 
apparent horizon). 

Penrose's heuristic argument for a future apparent horizon in this conjecture 
is described in more detail in [4j and [20] but roughly goes as follows: If, as is 
generally thought, asymptotically flat spacetimcs eventually settle down to a 
Kerr spacetime [12] , then in the distant future inequality [10] will be satisfied 
since explicit calculation verifies this fact for Kerr spacetimes, where A is the 
area of the event horizon. Given that some energy may radiate out to infinity, 
the total mass of these slices of Kerr may be less than the original total mass. 
Also, by the Hawking area theorem [10] (made more rigorous in [7|), and thus 
by the cosmic censor conjecture [24] as well, the area of the event horizon is 
nondecreasing in the spacetime evolution. Hence, this leads us to conjecture 
inequality [TO] in the initial Cauchy data slice, but where A is the total area 
of the event horizons of all of the black holes. The problem, though, is that, 
unlike apparent horizons, event horizons are not determined by local geometry 
but instead are defined in terms of which points in spacetime can eventually 
escape out to infinity along future directed time-like curves. Thus, in principle, 
there is no way to know which points this includes without looking at the entire 
evolution of the spacetime into the future. However, in [25] Penrose argued 
using the cosmic censor conjecture that future apparent horizons, which are 
defined in terms of local geometry, must be enclosed by event horizons. Thus, 
the area of E serves as a lower bound for the total area of the event horizons [16] , 
[13] ■ and the Penrose conjecture follows. This same argument, but run in the 
opposite time direction, yields the same conjecture for past apparent horizons 
as well. Thus, in the conjecture one could replace "future apparent horizon" 
with simply "apparent horizon." 

It is also important to note, as Penrose did originally, that a counterexample 
to the Penrose conjecture would be a very serious issue for general relativity since 
it would imply that some part of the above reasoning is false. The consensus 
among many is that the cosmic censor conjecture is the weakest link in the above 
argument. If the cosmic censor conjecture turns out to be false, and naked 
singularities (singularities not enclosed by the event horizons of black holes) 
do develop in generic spacetimes, then this would present a very interesting 
challenge to general relativity as a physical theory. 
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However, like the positive mass theorem, setting k = yields another beau- 
tiful geometric statement about manifolds with nonnegative scalar curvature, 
which is known to be true. 

Theorem 3 (The Riemannian Penrose Inequality, Bray, 2001 f^) 
Suppose that the Riemannian manifold {M^,g) is complete, has nonnegative 
scalar curvature, and is Schwarzs child at infinity with total mass m in a chosen 
end. If T,'^ € S is a zero mean curvature surface, then 

where A is the area of the outermost minimal area enclosure Y? = dU^ of 
TP'. Furthermore, equality occurs if and only if {M^ \ U^,g) is isometric to the 

Schwarzschild metric (ll^ \ B„/2(0), (l + ^j'' {dx^ + dy^ + dz^)j . 

In 1997, Huisken-Ilmancn proved a slightly weaker version of the above result 
with the modification that A is the area of the largest connected component of 
the outermost minimal area enclosure of and with the additional assumption 
that H2{M^) = 0. (This last topological condition can be replaced by assuming 
that TP is already a connected component of the outermost minimal area surface 
of {M^,g) by Meeks-Simon-Yau [22].) Their method of proof, first proposed 
by the theoretical physicists Geroch [S] and Jang-Wald ^^l, uses a parabolic 
technique called inverse mean curvature flow. Starting with a connected zero 
mean curvature surface, Huisken-Ilmancn found a weak definition of inverse 
mean curvature flow, where the surface is flowed out at each point in {M^,g) 
with speed equal to the reciprocal of the mean curvature of the surface at that 
point, for almost every surface in the flow. Then they showed that the Hawking 
mass of the surface is nondccrcasing under this flow, equals the right hand side 
of the Riemannian Penrose inequality initially, and limits to the left hand side 
of the Riemannian Penrose inequality as the surface flows out to large round 
spheres going to infinity. Both the physicists' insight into proposing this idea 
and the mathematicians' cleverness at generalizing the argument to something 
which could be made rigorous are remarkably beautiful. 

Bray's proof also involves a flow, but of the Riemannian manifold {M^,g). 
The flow of metrics stays inside the conformal class of the original metric and 
eventually flows to a Schwarzschild metric (shown as the case of equality metric) . 
The conformal flow of metrics is chosen so as to keep the area of the outermost 
minimal area enclosure of T constant. Also, the total mass of the Rieman- 
nian manifold is nonincrcasing by a clever argument (flrst used by Bunting and 
Masood-ul-Alam in [6]) using the positive mass theorem after a reflection of 
the manifold along a zero mean curvature surface and a conformal compacti- 
fication of one of the resulting two ends. Then since the Schwarzschild metric 
gives equality in inequality [Tl] the inequality follows for the original Riemannian 
manifold {M^,g). 
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All three systems of equations discussed in this paper which imply the Pen- 
rose conjecture are based on a new geometric identity which we call the gen- 
eralized Schoen-Yau identity. The identity is proved in section [SJ but with 
the lengthy computations relegated to the appendices for readability. This new 
idenity is a generalization of equation 2.25 in Schoen-Yau's paper [57]. The orig- 
inal Schoen-Yau identity was used to reduce the positive mass theorem to the 
Ricmannian positive mass theorem by solving a p.d.e. called the Jang equation. 
For all three systems, our technique will involve a generalization of the Jang 
equation to solving a system of two equations, the first of which is a generalized 
Jang equation in all three cases. Rather than spending time explaining the Jang 
equation, we will go straight to our proof since the Jang equation appears as a 
special case of our method (which for future reference is the case </> = !). 

As a final comment, the Penrose conjecture can be generalized to a statement 
about Cauchy data on n-manifolds motivated by considering (n-l-l)-dimensional 
spacetimes, where n > 3. In fact, the positive mass theorem was proved by 
Schoen-Yau in dimensions n < 7 and by Witten in any number of dimensions, 
but with the additional assumption that M" is spin. The Ricmannian Penrose 
inequality was proved by Bray [2] in dimension 3 using a proof which that author 
and Dan Lee [5] have generalized to manifolds in dimensions n < 7, and in a 
slightly weaker form by Huisken-Ilmanen |14| in dimension 3. Since we will be 
reducing the general case of the Penrose conjecture to the Ricmannian Penrose 
inequality, the techniques presented here have, at a minimum, the potential to 
address the Penrose conjecture for manifolds with dimensions n < 7. However, 
we will focus on 71 = 3 for simplicity. 

2 The Case of Equality 

In this section we carefully study the case of equality of the Penrose conjecture 
for the obvious reason that all of our estimates used to prove the conjecture must 
give equality in these cases. Also, we want to make sure that our techniques 
apply to all of the case of equality examples as a check that we are not making 
unjustified assumptions. 

We refer the reader to [23] for a discussion of the Schwarzschild spacetimc in 
Kruskal coordinates and follow those conventions (except for the names of the 
two functions a and /3 defined in a moment). Understanding the Schwarzschild 
spacetime in Kruskal coordinates is essential since this is the simplest global co- 
ordinate chart for the spacetime. In Kruskal coordinates, the entire Schwarzschild 
spacetime is expressed as the subset uv > —2m/e of x S*^ with coordinates 
{u,v,a S S*^) and line element 

2f3{r)dudv + r^da^, (12) 

where da^ is the standard round unit sphere metric on S'^, r > is a function 
oi u,v determined by 

uv = a{r) = {r- 2m)e(''/2™'-\ 
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and 

/3(r) = (8mVr)ei-(''/2™'. 

The first quadrant region described by u,v > is defined to be an exterior 
region and is isometrie to 



(^Rx(R^\B2™(0)),-(^l-— Jdt^+(^1- — ) dr^ + r^da") (13) 
under the isometry 



which we leave as an exercise for the interested reader to check. Note that 
we have now defined three different coordinate chart representations for the 
Schwarzschild spacetime, the two above in eauations ll2l and ll3[ and our original 
one in equation [2l 

A key point is that two of these coordinate chart representations of the 
exterior region of the Schwarzschild spacetime are written in the form of a 
static spacetime. For example, using the coordinates in equation ll31 the exterior 
region can be expressed as 

(R X A/'\ + g) , (14) 

where 



0^ = 1- 

which of course gives us 



2m 



2m 2m(j)'^ 
and r — 2m = 



1 - (/)2 1-02 

Hence, if we think of a sHce of the static spacetime expressed as the graph of 
t = f{x) in the static spacetime, a; S M^, we see that 

f = 2m\og{v/u), (15) 

W2°-P(^T3^j (16) 

The reason that these last two equations are important is that it allows 
us to understand the behavior of / and cj) as they approach the boundary of 
the exterior region {x \ u > 0,v > 0} oi the Schwarzschild spacetime. Our 
slice (intersected with the exterior region of the Schwarzschild spacetime) has a 
future apparent horizon boundary if it = everywhere on the boundary, a past 
apparent horizon boundary if w = everywhere on the boundary, and a future 
and past apparent horizon boundary if u = w = everywhere on the boundary. 

The mixed case where u = on part of the boundary and w = on the rest 
of the boundary does not represent a traditional apparent horizon boundary. 
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However, we note that whenever this boundary is area-outerminimizing, we are 
in fact in a case of equahty of the Penrose conjecture. This observation helps 
motivate the definition of a generaHzed apparent horizon in the next section. 

Also, while the u ~ level set and the v = level set on are both smooth 
(since the gradients of u and v on AI^ are never zero since is space- like), 
the boundary of {x € | u{x) > 0, v{x) > 0} in need not be smooth since 
the zero level sets of u and v do not need to intersect smoothly. However, when 
the boundary has corners it is never area outerminimizing and thus not a case 
of equality of the Penrose conjecture. 

We also note that apparent horizons outside of the exterior region, say 
with u = but with w < on part of the apparent horizon, are not area- 
outerminimizing since they have negative mean curvature at some points. Con- 
sequently, these last apparent horizons are enclosed by surfaces with less area 
and are therefore not cases of equality of the Penrose conjecture either. 

Kruskal coordinates reveals that the Schwarzschild spacetime is smooth on 
the boundary of the exterior region and certainly docs not have any singularities 
there. However, static coordinate representations of the Schwarzschild space- 
time have coordinate chart singularities there (which do not represent anything 
geometric or physical) . Hence, while Kruskal coordinates u and v are smooth on 
any slice, even up to the apparent horizon boundary, / and (f> are not necessarily. 

In fact, we see that / goes to ±oo logarithmically at the apparent horizon 
boundary typically (when u or v goes to zero and the other stays positive) . Also, 
(j)'^ vanishes on the apparent horizon boundary only linearly if either w or i; is 
strictly positive, which means that the derivative of (f) is going to oo. However, 
in the future and past apparent horizon boundary case where u, v both go to 
zero, then 0^ vanishes quadratically and tp is smooth up to the boundary. It 
is also true that / is smooth up to the boundary in this case by L'Hopital's 
rule since u and v, which equal zero on the future and past apparent horizon, 
have nonzero derivatives there (since the hypersurface is space-like). These 
observations are helpful since we will be dealing with slices of the exterior region 
of the Schwarzschild spacetime viewed in static coordinates for the rest of this 
paper. 

3 Generalized Apparent Horizons and the Gen- 
eralized Penrose Conjecture 

In this section wc describe the most general version of the Penrose conjecture 
that we believe to be true. Naturally this is an important question to consider 
since proofs of a conjecture may be more easily found when the most natural 
version of the conjecture is understood. 

Definition 7 Define the smooth surface Yp' £ S in (M^, g, k) to he a generalized 
apparent horizon if 

i?s = |trs(fc)| (17) 
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and a generalized trapped surface if 

< |trs(fc)|. 



(18) 



In terms of the mean curvature vector H of in the spacetime, a generahzed 
trapped surface is one where H is not strictly inward space-Hke anywhere on 
S. Also note that this definition of a generahzed apparent horizon does not 
need a globaUy defined future directed unit normal to since the definition 
is unaffected by a change of sign of the second fundamental form k. A related 
class of surfaces, referred to as "^-surfaces" , appears in a different context in 

Referring to the previous section, note that any smooth slice of the 
Schwarzschild spacetime which smoothly intersects (which is often not the case) 
with the boundary of the first quadrant {m > , w > 0} of Kruskal coordinates 
intersects in a generalized apparent horizon. These generalized apparent hori- 
zons also give equality in the Penrose conjecture, so it is natural to include them 
in the statement of a generalized Penrose conjecture. Also note that traditional 
apparent horizons, if they are not already generalized apparent horizons, are at 
least always generalized trapped surfaces. 

Another consideration which leads to this definition of generalized apparent 
horizons is the case when a surface with multiple connected components is a 
future apparent horizon on some connected components and a past apparent 
horizon on the others. While Penrose's original heuristic argument does not 
apply to this surface, the techniques that we develop in this paper seem to 
apply perfectly well. Thus, we would like a generalized Penrose conjecture 
which includes this case as well. 

After a talk on generalized apparent horizons by the first author at the 
Niels Bohr International Academy's program "Mathematical Aspects of General 
Relativity" in April 2008, Robert Wald posed the following insightful question: 
Is it possible for generalized trapped surfaces to exist as boundaries of space-like 
slices of Minkowski space? (A similar query was posed by Mars and Senovilla 
in |21|.) This question raises the issue of whether or not generalized trapped 
surfaces always yield a positive contribution to the ADM mass, which of course 
is a prerequisite for a generalized version of the Penrose conjecture. If one could 
find a generalized trapped surface which was the boundary of a space-like slice 
of Minkowski space, then the total mass of the slice would be zero, making a 
Penrose-type inequality for the surface impossible. 

In response to this question, the second author of this paper showed that 
no such generalized trapped surface in Minkowski space exists }18) . Further- 
more, he showed that Witten's proof of the positive mass theorem also works 
for asymptotically flat manifolds with generalized trapped surface boundary and 
gives a positive lower bound on the total mass. This result suggests that gen- 
eralized trapped surfaces and generalized apparent horizons have some physical 
significance in that such surfaces, along with nonnegative energy density /x > | J| 
everywhere in the spacetime, always imply that the total mass is positive. Find- 
ing the best possible lower bound on the total mass motivates conjecturing a 
generalized Penrose inequality. 
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In addition, a discussion between the first author and Tom Ilmanen led to two 
more conjectures about generahzed apparent horizons and generaUzed trapped 
surfaces, which are known to be true in the special case fc = 0. We are pleased 
that Michael Eichmair [Hj has announced proofs of these two conjectures (except 
for the topological part of conjecture [3]) using elliptic techniques (whereas Ilma- 
nen's original ideas used parabolic techniques). We omit the n = 2 case in these 
next two conjectures because they are less relevant for our present purposes, 
but we understand that Eichmair's results apply there as well. 

Conjecture 2 (Tom Ilmanen, 2006) 

Given complete, asymptotically flat Cauchy data {M^^,g,k), 3 < n < 7 , with a 
generalized trapped surface then there exists a unique outermost general- 

ized trapped surface E which is a generalized apparent horizon. 

Conjecture 3 (Tom Ilmanen, 2006) 

Furthermore, S is strictly area outerminimizing ( every other surface which en- 
closes it has larger area), and for n = 3, the region exterior to S is diffeomorphic 
to R'^ minus a finite number of disjoint closed balls. 

The above conjecture is a generalization of Mccks-Simon-Yau [22], which is 
the case when k = 0. The topological conclusions of this last conjecture, like the 
original Meeks-Simon-Yau result, make this conjecture particularly interesting 
for its own sake as well as important for the Jang-IMCF system of equations 
we will describe later in the paper. We also encourage the reader to study the 
related theorems of Andersson and Metzger [1] on future and past apparent 
horizons, which are relevant for this discussion. 

All together, these considerations lead us to make the following generalized 
Penrose conjecture. Since traditional apparent horizons are always generalized 
trapped surfaces, this conjecture implies the original Penrose conjecture. 

Conjecture 4 (The Generalized Penrose Conjecture) 

Suppose that the Cauchy data {M^,g,k) is complete, satisfies the nonnegative 
energy density condition fi > \J\, and is Schwarzschild at infinity with total 
mass m in a chosen end. IfTF'&Sisa generalized trapped surface, then 

where A is the area of the outermost minimal area enclosure T? = dU^ of Y? . 
Furthermore, equality occurs if and only if {M^ \U^,g,k) is the pullback of the 
Cauchy data induced on the image of a space-like embedding of M^\U^ into the 
exterior region of a Schwarzschild spacetime (which maps to a generalized 
apparent horizon). 

We note that this conjecture is true when fc ~ by [5]. In this case, S has 
nonpositive mean curvature and acts as a barrier to imply the existence of an 
outermost minimal area enclosure of E which is minimal. 
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It is important to note that conjectures [Hand [3] imply that the generaUzed 
Penrose conjecture (and hence the original Penrose conjecture) follows from the 
following important case of the generalized Penrose conjecture. 

Conjecture 5 (The Generalized Penrose Conjecture - Outermost Case) 
Suppose that the Cauchy data {M'^,g,k) is complete, satisfies the nonnegative 
energy density condition fi > \J\, and is Schwarzschild at infinity with total 
mass m in a chosen end. Suppose also that = dU^ € S is a strictly area 
outerminimizing generalized apparent horizon, that no other generalized trapped 
surfaces enclose it, and that \ is diffeomorphic to R'^ minus a finite 
number of disjoint closed balls. Then 



where A is the area of . Furthermore, equality occurs if and only if (M^ \ 
U^,g,k) is the pullback of the Cauchy data induced on the image of a space- 
like embedding of \ into the exterior region of a Schwarzschild spacetime 
(which maps to a generalized apparent horizon). 

Most of the remainder of this paper will focus on proving the above con- 
jecture. Naturally, when attempting a difficult conjecture, it makes sense to 
consider the simplest case which still captures the essential subtleties of the prob- 
lem, and the above conjecture is arguably that case. We will focus on the three 
dimensional case in this paper, but the above conjecture is the same in higher 
dimensions up to and including seven, but without any hypothesis on the topol- 
ogy of M" \ U", where the conjectured inequality becomes m > c„^("~^^/("~^) 



Furthermore, the condition of not having any generalized trapped surfaces 
outside of S may turn out to be very important. There is reason to believe 
that the generalized Jang equation, which is a p.d.e. we will define later in the 
paper which is central to all of our approaches, may blow up on surfaces with 
\Hy:\ = |trs(fc)| (which are clearly generalized trapped surfaces). Since we are 
leaving the existence theory of our proposed p.d.e. 's which imply the Penrose 
conjecture open, those who arc considering studying these existence theories 
will have to understand this possible behavior carefully. 

4 Proof of the Penrose Conjecture in a Special 
Case 

In this section we will prove the Penrose conjecture, conjecture [U with two 
extra assumptions, and show how the conjecture follows from the Ricmannian 
Penrose inequality, theorem [3l This special case, where a correct approach is 
quite clear, will help us motivate the general case which is not so obvious. 

A major hint in the statement of the Penrose conjecture is the case of equal- 
ity. Since the Penrose conjecture is an equality for any slice (space-like hyper- 
surface) of the exterior region of the Schwarzschild spacetime with an apparent 




(20) 



see [5]. 
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horizon boundary, we know that all of our techniques must preserve this equality 
in every estimate we derive. 

On the other hand, if we are given some Cauchy data {M^, g, k) which comes 
from a slice of the Schwarzschild spacetime, it may be difficult to recognize it 
as such. However, our techniques must absolutely be able to recognize these 
Cauchy data as the instances where we get equality in all of our inequalities. 

More generally, suppose (M^,g, k) comes from a slice of the static spacetime 

(R X A/3, ~(t>'^df + g) , (21) 

where is a real- valued function on M and g is some other Ricmannian metric 
on M . Notice that the Schwarzschild spacetime can be expressed in the form of 
equation 1211 However, while the Schwarzschild spacetime is vacuum (meaning 
it has zero Einstein curvature and consequently zero Ricci curvature), we are 
making no such requirement on (R x Af'^, —(fP'dP + g). 

Given a real- valued function f on M , define the graph map 

F : M ^Rx M (22) 

where F{x) = {f{x),x). Then a short calculation reveals that the puUback of 
the induced metric on the image of _F in a coordinate chart is gij — 4>'^ fifj, so 
setting 

guarantees that the puUback of the induced metric on the image of the graph 
map F is precisely g. A similar type of calculation (but which is much longer 
and so is carried out in the appendices) yields that the puUback of the second 
fundamental form of the image of the graph map F in the static spacetime to 
(M3,5) is 

, _ (t>Hessijf + (t)ifj + fi(l)j 
y ^ 7772 ' ^ ' 

where subscripts on / and (j) represent coordinate chart partial derivatives and 
the Hessian of / is taken with respect to the metric g (or the Levi-Civita con- 
nection of g if one prefers) . These considerations lead us to the following special 
case of the Penrose conjecture which has an elegant and relatively short proof 
using the Gauss-Codazzi identities and the Ricmannian Penrose inequality. 

Theorem 4 The Penrose conjecture as stated in conjecture [7] follows for 
{M^,g, k) if there exist two smooth functions f and (j) on such that 



^ 02|wf|2^l/2 

and 



fc.,=fe.,^ ^''^^^'^^ + '^'^+^^^ outs^deofl: (24) 



= on S, (25) 
where (f) > outside of S and f has compact support. 
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Proof: We will reduce the Penrose conjecture on {M^,g,k) to the Rieman- 
nian Penrose inequality on {M^,g). To do this we need to show that 

• the scalar curvature ^ of 5 is nonncgativc and that 
o S has zero mean curvature H in (M^,g). 

Then the fact that g measures areas to be at least as large as g docs implies 
that the area of any surface in {M^,g) is at least as large as the area of that 
same surface in {M^,g). Thus, 



where Sg and Sg are the outermost minimal area enclosures of S in (M'^,g) 
and {M^,g), respectively. Since / has compact support, the masses of the two 
manifolds are the same. Then by the Ricmannian Penrose inequality 



which proves that Penrose conjecture on {M^,g,k). Thus, all that is left to 
prove are the two bullet points (•) and (o). 

Proof of (o): Since cjij = gij+<t>^ fifj and = on S and (j) and / are smooth, 
the two metrics are the same up to first order on E. But the mean curvature 
of a surface, which is the main term in the first variation of area formula, only 
depends on the metric and the first derivatives of the metric. Hence, H ~ H . 

Since S is an apparent horizon, H = ±trs(A;). But since </> = on S, 
derivatives along S of are zero as well, so our assumption on the special form 
of k in equation [Ml implies that trs(fc) = 0. Hence, 



Proof of (•): Working inside of the static spacetime in equation [21] let n be 
the future pointing normal vector to the image of under the graph map F 
from equation [55] and let n be the future pointing normal vector to M'^ viewed 
as the t = Q slice of the static spacetime. Then these two vector fields on 
hypersurfaces can be extended to the entire spacetime by requiring that these 
extended vector fields are invariant under translation in the time coordinate 
(which is an isometry of the spacetime). 

The trick is to compute G{n, n) using the Gauss-Codazzi identities, but in 
two different ways. We are given the nonnegative energy density condition on 
{M^,g,k) that /i > |J|. Since we are in the very special case that k actually 
equals the second fundamental form h of the graph, {M^,g,h) has /i > |J| too. 
This is equivalent to saying that G{n, w) > for all future time-like vectors w 
in the spacetime. Letting w = n thus implies that 




H = H = trE(fc) 0. 



G{n,n) > 0. 



(26) 
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On the other hand, applying the Gauss-Codazzi identities to the t = shce of 
the static spacetime gives us 



{8n) p ^ G{n,n) = (i? + tr(p)2 - ||pf )/2 
(Stt) J = G{n,-) = div (p - tr(p).g) 

where p is the second fundamental form of the t = slice, which of course is 
zero by the time symmetry of the spacetime. Hence, (Stt)/! = R/2 and J = 0. 
Thus, if we let 

n = an + (vector tangent to t = slice), 
where a is a positive function on M, we have that 

G{n, n) = aG{n, n) = aR/2. (27) 

But G is symmetric, so by inequality 1261 -R > 0, which completes the proof of 
(•) and the proof of the Penrose inequality in this special case. 

The case of equality of the above theorem would follow from conjecture [7] 
in section [T) We refer the reader to that section for discussion on the case of 
equality since the main purpose of this section was to motivate the identities 
computed in the next section. 



5 The Generalized Schoen-Yau Identity 

The proof of the Penrose conjecture in the special case presented above suggests 
how the Gauss-Codazzi identities can be used to compute a formula for the scalar 
curvature Roig = g~\- (f^dp in terms of the scalar curvature R of g, the graph 
function /, and the warping factor 0. In this section we will derive this formula 
and then show how this formula leads to an identity central to our approach to 
the Penrose conjecture. 

From this point on we will abuse terminology slightly and always refer to 
the image of the graph map F{M) simply as M and the t = Q slice of the 
constructed spacetime as M . This notation is convenient since then (M, g) and 
(M, g) are space-like hypersurfaces of the spacetime (R x M^, —(fidt^ -\- g) . Let 
TT : M M be the projection map TT{f{x),x) = (0,x) to the t = slice of the 
spacetime. 

Establishing some notation, let Bq = dt and {di} be coordinate vectors 
tangent to M. Define 

^^ - a, + m (28) 

to be the corresponding coordinate vectors tangent to M so that tt, (9^) = di. 
Then in this coordinate chart, we have that 

9ij =913 ~ 4>^fif3- (29) 

It is convenient to write 

g'^ =g'j +v'v^ , (30) 
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where 

" (l-02|d/|2)l/2 (l+^2|rfJ|2)l/2- ^ ) 

We also define 

v = v'^^ and w = w^^i (32) 
so that TT, (w) = v, and observe the useful identity 

{l^ct^^\df\l)-{l + ^^\df\l)^l, (33) 

which is evident by looking at the ratios of the volume forms. See appendix [Cl 
for more discussion on these calculations. 

In this paper we use the convention that a barred index (as in p above) 
denotes an index raised (or lowered) by g as opposed to g. That is, /* = 
g^-' fj, where as usual fj = df /dxj in the coordinate chart. In general, barred 
quantities will be associated with the t = slice {M, g) and unbarred quantities 
will be associated with the graph slice {M,g). 

In appendix [P] we compute that the second fundamental form of the graph 
slice (M, g) in our constructed static spacetime is 



, _ '/>Hess,j-/ + {fi(l)j + (/),;/j) - d(l))gf,fj 

- (l-02|d/||)V2 ^-^^^ 

_ 0HesSjj/ + (/i^j + ^Jj) ^^^^ 



(l+02|d/|2)l/2 
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which we Ust now for future reference. 

Finally, we extend h and k trivially in our constructed static spacetime 
so that h{dt, ) = = k{dt,-) and such that these extended 2-tensors equal 
the original 2-tensors when restricted to M. Note that this gives h{di,dj) = 
h{di,dj), so we can call this term hij without ambiguity. The same is true for kij 
and components of 1-forms like /, and (pi. However, we remind the reader that 
the Hessian of a function, which is the covariant derivative of the differential of 
a function, depends on the connection and hence the metric since we will always 
be using the respective Levi-Civita connections on {M,g) and {M,g). 

Now we are ready to proceed to compute a formula for R. It is a short 
calculation to verify that, in the constructed spacetime, 

(n,n) = -(1 - <l^'\df\lr'/' = -(1 + d^'ldfll)'/' (36) 

Thus, 

n = (1 + <f>'\df\l)'/^n + tang^ap^(n) (37) 
where another short calculation reveals that 

tangraph(") = -^f'^j - -0V/. (38) 

As in the previous section, the trick is to compute G(n, n) two different ways 
using the Gauss-Codazzi identities. As before, applying these identities to the 
t = slice (M, g) of the constructed spacetime gives us 

G{n,n) = (l + 02|rf/|2)i/2G(^^^) 

= (l + 02|^^|2)l/2.^/2 

since the t = slice has zero second fundamental form. On the other hand, 
applying the Gauss-Codazzi identities to the graph slice {M, g) yields 

G(n,n) = (l + 02|d/|2)i/2G(n,n) + G(n,tang^api^(n)) 

= {l + ^^\df\l)^'^[R+{ir,hf - \\h\\l]/2 
-fdiv(/i-(tr<,%)(-0V/). 

Combining the two previous equations, we get our first desired result 

R = R+ {tvghf - \\h\\l + 2{d{tvgh) - 6:iY{h)){v). (39) 

Of course, what we are given in the hypotheses of the Penrose conjecture is 
that /i > \ J\g, where 

(8^) = G(n,n) = (i? + trg(A:)2 - ||fcf)/2 
(Stt) J = G(7i, •) = div(A;) - d(trg(fc)), 

for some symmetric 2-tensor k. Hence, 

R = l%7:{^,^J{v)) + {trghf-{trgkY-\\h\\1 + \\k\\l 

+2v{ti-gh) ~2v{tigk) -2div{h){v) + 2div{k){v). (40) 
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Note that fi — J{v) > since \v\g < 1. Hence, as we saw in the previous 
section, if we can choose a </> and an / so that h = k, then we immediately 
get that ^ > 0. However, we are interested in investigating if a more general 
relationship between h and k can give a similar result. 

Our procedure is to convert our formula for R to an expression in terms of the 
g metric. Arguably g is more natural than g since it is the metric induced on the 
t = slice of the static spacetime. To perform the conversion, we need several 
identities for arbitrary symmetric 2-tensors k which are proven in appendix |D] 
and which we list here. 



Identity 1 



trg{k)f - \\k\\l = {trgkf - \\k\\l + 2k{v,v)tr,k - 2\k{v, 



Identity 2 
Identity 3 
Identity 4 



v{trgk) — v{trgk + k{v, v)) 



dm{k){v) = dw(fc)(w) + (Vsfc)(w,w) -2|w|^ fc I w,^ 



+ {h{v, •), k{v, ■))g + 2h{v, v)k{v, v) + {trgh)k{v, v) 



h,j + vjh[v, ■)j -f- 



Identity 5 



Identity 6 

div{k){v) = div{k{v, •)) — {h, k)g — {h{v, •), k{v, ■))g + k ^u, 
Identity 7 

(Vc/c)(w, v) = v{k{v, v)) - 2{h{v, ■), k{v, ■))g - 2h{v, v)k{v, v) + 2\v\l k (^v, ^ 
Identity 8 



div(k){v) = div{k{v,-)) + v{k{v,v)) + k yv,- 

-{h, k)-g - 2{h{v, •), k{v, ■))-g + {tr-gh)k{v, v) 

Identities [T] and are short calculations. Identity [3] is used in the proof of 
identity ID Plugging identities [6] and [7] (which are proved using identity [5]) into 
identity [4] results in identity [H Finally, plugging identities [U [2l and [8] into our 
formula for R results in the main identity of this paper. 
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Identity 9 (The Generalized Schoen-Yau Identity) 

R = 167r{fi~Jiv)) + \\h-k\\l + 2\q\l~p^i(f>q) 

+ {trghY - [trgkY + 2v{trgh - trgk) + 2k{v, v){trgh - tVgk) 

where 

q = h{v, •) - k{v, ■) = h{v, ■) - k{v, •) . 

Note that the two definitions of q exist on the entire constructed static spacetime 
and are equal since both h and k are extended trivially in the constructed static 
spacetime. We also observe that 

\div{(j)q) ^ divsT(g), 

where div^T- is the divergence operator in the constructed static spacetime. 

In the special case that </>=!, the above identity was derived by a different 
method by Schoen-Yau as equation 2.25 of [57j (in fact the procedure in [57] 
may also be used to obtain identity [S] and will be presented in a future paper) . 
In that paper, they used the Jang equation, 

= tTg{h - k) 

to reduce the positive mass theorem to the Ricmannian positive mass theorem. 
While imposing the Jang equation in the special case that 0—1 does not imply 
that _R > as would be most desirable, R > 2\q\^ — 2dYv{q) implies that there 
exists a conformal factor on g such that the conformal metric has nonnegative 
scalar curvature and total mass less than or equal to that of g and g. Then the 
Riemannian positive mass theorem applied to the metric conformal to g implies 
the positive mass theorem on (M, g). This approach docs not quite work for the 
Penrose conjecture because the conformal factor needed to achieve nonnegative 
scalar curvature changes the area of the horizon in a way which is difficult to 
control. 

6 The Generalized Jang Equation 

All of the approaches to the Penrose conjecture that we consider in this paper 
use the generalized Schoen-Yau identity. This identity plays a central role in the 
remainder of our discussions because it directly relates the nonnegative energy 
condition on {M^,g,k) (which implies that fj, > J{v) since \v\g < 1) to the 
scalar curvature of {M^,g). 

Furthermore, this generalized Schoen-Yau identity strongly motivates the 
generalized Jang equation, 

= trg{h~k), (41) 
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which on the original manifold (M^, g) with Cauchy data {M^,g, k) is the equa- 
tion 

c^^fP \ ( <t>Hess,,,f + c^,fj + 



l + C^'\df\lJ V (l + 02|d/|2)V2 



(42) 



when one substitutes the formulas for h and g'^ in a coordinate chart. (In this 
paper we adopt Einstein's convention that whenever there are both raised and 
lowered indices, summation is implied, so the above formula is a summation 
over i,j both ranging from 1 to 3.) 

Of course the original Jang equation, which again is the special case (j){x) = 1 , 
only had one free function, /, whereas the generalized Jang equation has two 
free functions, / and (j). Hence, to get a determined system of equations, we 
need to specify one more equation. Later in the paper we will propose various 
choices for this second equation, but our choice for the first equation will always 
be the generalized Jang equation above. 

Once the generalized Jang equation is specified, the generalized Schoen-Yau 
identity simplifies greatly to 

R = 16ni^i - J{v)) + \\h - k\\l + 2\q\l - |m09). (43) 

It is important to note that the first three terms of the right hand side of the 
above equation are all nonnegative since /i > \J\g and \v\g < 1. 



6.1 Boundary Conditions 

Examining the case of equality slices of the Schwarzschild spacetime described 
in section[2] leads us to propose the following boundary conditions on generalized 
apparent horizons. At a minimum, these boundary conditions are satisfied al- 
most everywhere for slices of the exterior region of the Schwarzschild spacetime 
with generalized apparent horizon boundaries. 

Boundary Conditions on Generalized Apparent Horizons 

Given a generalized apparent horizon E with mean curvature H^: = | tr|,(fc)| and 
outward unit normal v in {M^,g), we require that (j) = and 

{i^,v)g = sign{trl,ik)) (44) 

on S, where as usual 

" (l + 02|rfJ|2)l/2 

and V is extended to the boundary E by continuity. 

Note that these boundary conditions are consistent with / blowing up to +00 
where tr|,(A:) < 0, blowing down to —00 where tr|;(/c) > 0, and / staying 
bounded where tr|,(A:) = on E. 
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The hope is that these or similar boundary conditions imply that S, which 
was a generalized apparent horizon in {M^,g), becomes a minimal surface with 

i?s = 

in {M^,g). We discuss the general case of this question in appendix lEl For 
now, we observe two important special cases. 

The first important special case is when S is a traditional apparent horizon, 
either future or past, and Hy: > 0. If we also assume that / goes to ±00 on 
each connected component of E in a reasonable fashion, then the level sets of / 
converge to E. The formula for the mean curvature of the level sets of / in the 
new metric g is 

H = (i+(/)2|d/n-i/2ij 

since g ~ g + (f^dp docs not change the metric on the level sets of /. stretches 
lengths perpendicular to the level sets of / by a factor of (1 + and 
by the first variation formula for area. Then if we assume that / and (p behave 
similarly to the case of equality slices of Schwarzschild, we get the following 
lemma. 

Lemma 1 Suppose that {M^,g,k) has a smooth interior boundary E which 
is a future [past] apparent horizon with Hy, > 0. Then if f blows up [blows 
down] logarithmically, \df\ blows up asymptotic to 1/s, and 0^ goes to zero 
asymptotic to s (where s is the distance to E in {M^,g)), then the limit of the 
mean curvatures H of the level sets of f in {M^,g) is zero. 

The second important special case is the case where the boundary is a future 
and past apparent horizon. In this case, based on the case of equality slices of 
Schwarzschild, we expect / to stay bounded and smooth and (p to stay smooth 
as well. 

Lemma 2 Suppose that {M^,g,k) has a smooth interior boundary E which 
is a future and past apparent horizon (which by definition has H-^ — and 
tr|,(fc) = 0). Then if f is bounded and smooth and (f> is smooth and equals zero 
on E, then H^: = i^s = 0. 

The proof of this lemma appeared in this paper already in section ID The point 
is that since g = g + (p'^df^, both metrics g and g are the same up to first order 
on E since / and ^ are smooth and (/> = on E. Then since the mean curvature 
is only a function of the metric and first derivatives of the metric, the two mean 
curvatures are equal, and since i7s = 0, both are zero. 

6.2 Blowups, Blowdowns, and Outermost Horizons 

One phenomenon of the original Jang equation ((/) = 1) is that / can blowup to 
00 on future apparent horizons or blowdown to —00 on past apparent horizons. 
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and this feature is still present in the generalized Jang equation (given plausible 
assumptions about the behavior of / and ^). More importantly, according to 
|27| , blowups and blowdowns of / with the original Jang equation can only occur 
on apparent horizons. 

An important question, then, is to understand when blowups can occur 
with the generalized Jang equation. Certainly blowups and blowdowns can still 
occur on traditional apparent horizons. A reasonable conjecture is that the 
blowup properties of the generalized Jang equation are the same as the original 
Jang equation as long as (j) is smooth and strictly positive. However, if (f> is 
allowed to go to zero, then we have already seen that / can have a mixture 
of blowup, blowdown, and bounded behavior on generalized apparent horizons 
in case of equality slices of the Schwarzschild spacetime. Those who study the 
existence theories of the systems of equations proposed in this paper will need 
to understand these issues. 

A reasonable hope, however, is that as long as our boundary E is already an 
outermost generalized apparent horizon, so that no other generalized trapped 
surfaces enclose E, then / stays bounded away from S as long as (j) stays strictly 
positive. A relevant calculation which is useful for studying the question of when 
/ can blowup or blowdown is the following. 

Recall the standard identity 

A/ = Hess/(i., //) + i/s • Kf) + As/ 

for the Laplacian of a function in terms of the Laplacian of that function re- 
stricted to a hypersurface with mean curvature H and outward unit normal v. 
If we let S be any level set of /, then we get that 

tr|(Hess/) = T\df\gH^ 

for blowup and blowdown respectively. Thus, the generalized Jang equation 
implies that 

= trg{h~k) 

{h - k){v, v) 



tr|(/i-A:) + 



a 



4> tr|(Hcss/) _ , 5 /, ^ ih-k){v, v) 

(l + 02|d/|2)l/2 si J+ l + ^2|rfj|2 

-r '^l^/lg - tr^ (k) + 



l + 02|d/| 
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on level sets of /, where we have used the facts that v and v are collinear, 



Ml = 1 - 



l + 02|rf/|2' 

and the formulas for g'^ and hij from section [5] 

Lemma 3 When f is blowing up or blowing down, the term h{v, u) is bounded 
if(j)'^df is assumed to be smooth and nonzero in the limit up to the boundary and 
(j) = on the boundary, which is true in case of equality slices of Schwarzschild. 

Proof: Referrmg back to section [21 in a smooth slice of Schwarzschild (fp'df 
can be expressed in terms of the smooth Kruskal coordinate variables u, v. Since 
by equation [T5l 

and by equation [TBI 

(f)^ = uvj{uv) 

for some smooth function 7 7^ for 0^ < 1, we have that 

(ji^df = 2m'f{uv) {udv — vdu) . 

The fact that our case of equality slices of Schwarzschild are spacelike implies 
that \du\g ^ 7^ \dv\g, and the fact that we are assuming blowup or blowdown 
implies that exactly one of u, v is going to zero on the boundary. Hence, not 
only is ip^df smooth up to the boundary, it is also nonzero in the limit up to 
the boundary. 

Then the fact that h{v, v) is bounded up to the boundary assuming this 
smoothness follows from the short calculation that 

V.{cl>^df){v) 

which completes the proof of the lemma. 

Thus, referring back to our calculation before the lemma, given blowup or 
blowdown with behavior on / and <f> as seen in the case of equality slices of 
Schwarzschild in scction[21 (t>^ goes to zero linearly, \df\g goes to infinity like 1/s 
so that (fp'\df\'^g goes to infinity like 1/s, and h{v, v) stays bounded. Then since k 
is given to be smooth and therefore bounded, we conclude that the generalized 
Jang equation implies that 

= Ti?E-tr|(fc) 

on surfaces with this type of blowup or blowdown of /, which of course are the 
equations for future and past apparent horizons, respectively. 
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7 The Jang - Zero Divergence Equations 



Looking at equation l43l the most direct way to get > is to set div{(j>q) ~ 0. 
We will call the resulting system of equations, equations |46] and [47l the Jang 
- zero divergence equations. The following existence conjecture for these equa- 
tions implies the outermost case of the generalized Penrose conjecture, conjec- 
ture [HI using the Riemannian Penrose inequality, and is therefore an important 
open problem. 

Conjecture 6 Given asymptotically flat Cauchy data {AI,g,k) with an outer- 
most generalized apparent horizon boundary S, there exists a solution (/, ip) to 
the system of equations 

= trg{h~k) (46) 
= (47) 

with lim,_oo/(.T) = 0, 02|v/|2 = o(r-i), V((/.2 1 V/|2) = o{r-^), 
and lim^^^oo 0(a;) = 1, where g = g + 4>'^df'^ , 

(j)Hessf + {df ®d(p + d(j)® df) 



{l + c^^\df\lY/^ 



(48) 



q = h(v, •) — k{v, ■), and v = (pVf /{I + such that S has zero mean 

curvature in the g metric. 

The boundary conditions on / which lead to E having zero mean curvature in 
(M^ , g) are discussed in the previous section and in appendix [E] 

Also, we comment that while equation |47] is third order in /, subtracting 
derivatives of equation |46] can remove the third order terms of / in favor of Ricci 
curvature terms. While the resulting system has quadratic second order terms 
in / in the second equation, the system is degenerate elliptic. 

The above system may also be reduced to a system of 1st order equations 
by introducing new variables. If we let a = df, then the above system has 
a solution whenever the first order system with variables </) (a 0-form), a (a 
1-form), and /3 (a 2-form) 

^ da (49) 
= d/3 (50) 
= tig{h-k) (51) 
(j)q = d*/3 (52) 

has a solution, where d* is the d star operator with respect to the g metric 
which sends 2-forms to 1-forms. In these variables, g = g + (j^^ct^, 

_ (t>Va + {a ® d4> + dcf) ® a) . 

" (l + 02|^|2)l/2 ' y^-^) 

q = h{v, ■) — k{v, •), and v = 0(3/(1 -I- (j)'^\a\'^y^^ , where a is the dual vector to 
a with respect to g. 
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Theorem 5 Conjectures\^ and^ (defined below to handle the case of equality) 
imply conjecture\Bi the outermost case of the generalized Penrose conjecture. 



Proof: As was discussed in the previous section, the point of requiring E to be 
outermost in the above conjecture is so that / does not blowup or blowdown on 
the interior of A/. Then the method of proof assuming conjecture [6] is basically 
the same as the proof of the Penrose conjecture in the special case in section 
[31 The total mass of {M^,g) is the same as the total mass of {M^,g) since the 
total mass is defined in terms of the l/r rate of decay of the metrics which are 
equal since ||g — gUg = Also, since g measures lengths, areas, etc. to 

be greater than or equal to that measured by g, 



Hence, the Penrose conjecture on {M^,g, k) follows from the Riemannian Pen- 
rose inequality on {M^ ,g). 

Thus, all that remains is to show that the Riemannian Penrose inequality can 
be applied to (M'^, g). The existence theorem already gives us that H = Q,&o the 
last thing to check is that R > 0, which follows directly from the generalized 
Schoen-Yau identity and equation This proves the inequality part of the 
Penrose conjecture on {M^,g,k). 

In the case of equality of the Penrose conjecture, clearly we must have equal- 
ity in all of our inequalities. Since the case of equality of the Riemannian Pen- 
rose inequality is solely when g is the Schwarzschild metric which has zero scalar 
curvature R, equation 1431 gives us 



Since each of these three terms is nonnegative, each must be zero. Hence, 
k = h. If we could argue that c/> = 0o, where (j)o is the warping factor from 
the Schwarzschild spacetime, then we would have that k = h is the second 
fundamental form and g = g ~ (ji^df"^ is the induced metric of a slice of a 
Schwarzschild spacetime, as desired. 

However, there is a delicate point here. In fact, (j) does not have to equal 
(po for {M^,g,k) to be the Cauchy data from a slice of a Schwarzschild space- 
time. If = c(f>o for some constant c > 0, then defining dfo = cdf (which can 
be integrated to recover /o) implies that (/o,^o) and (/,</)) produce the same 
metrics and second fundamental forms. This may seem like a minor point at 
first, but in fact this statement is still true if (and only if which we leave as an 
exercise) dc = on the open region D where df ^ 0. Thus, c may be different 
constants on each connected component of D. Again, dfo, which is still closed, 
may be integrated to recover /p since the t = slice of Schwarzschild is simply 
connected. Thus, we have the following lemma. 

Lemma 4 // = ctpo, where dc ^ on {x \ df ^ 0}, then (M^,g,k) comes 
from a slice of the Schwarzschild spacetime. 



(54) 



0=l&^{^L-J{v)) + \\h-k\\l + 2\q\l 



(55) 
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To prove the case of equality of the Penrose conjecture then, we need to 
prove the hypotheses of the above lemma. Looking back at equation [55l we see 
that we must also have 



= ^ - J{v) ^ - I J|) + I J|(l - \v\) + {\J\\v\ - J{v)), 

where all norms arc with respect to g. Again, since each of the three grouped 
terms is nonnegative, all must be zero. Since \v\ < 1, the second term equaling 
zero implies that \J\ = so that the first term equalling zero implies that /i = 0. 
In the appendices we compute that in the static spacctime g — (j)'^dt^, 



(l-02|d/|2)l/2(,-, + ^V/) 



and that if i? = 0, 



J = G(n,-) 



Ric 



Hess(/) A(/) , 



(0V/, ■ 



where • is a tangent vector to the graph slice {M^ ^ g) in the first instance and its 
component tangent to the i = slice in the second. Then since the Schwarzschild 
spacetime has G = 0, = 0, A(/)o = 0, and 



Hess^o 

RlC = ; , 



J = in the case of equality implies the ovcrdctcrmincd equation (when df ^ 0) 
for (f) that 

'A0 A</.o~ 







Hcss0o Hcssf/) 



(V/,-). 



(56) 



Conjecture 7 Eauation \56] imvlies the hypotheses of lemma^ 



Clearly the hypotheses of lemma 2] imply equation 1561 but we need the 
converse to be true as well. Assuming conjecture IHl is true, a proof of conjecture 
[7] would finish the case of equality part of the outermost case of the generalized 
Penrose conjecture. 



8 Einstein-Hilbert Action Methods 

Equation [33] is a remarkable equation which deserves very careful consideration. 
Since we need a lower bound on the scalar curvature R of {M^,g), the only 
troublesome term in that equation is the last one, the divergence term. In the 
previous section, we dealt with this last term by setting it equal to zero. In 
this section, we make the natural observation that divergence terms can also be 
dealt with by integrating them. 
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Theorem 6 7/ g = g + (p'^df'^ on with boundary and Cauchy data 
{M^,g,k) satisfying the nonnegative energy condition fi > \J\, the generalized 
Jang equation 

= tTg(h - k) 

is satisfied, and f and <j) behave at infinity and on the boundary E such that 
eauations \60\ and \61\ are satisfied as expected, then 

R4)dV>0, (57) 

M 

where R is the scalar curvature of g and dV is the volume form of g. 

In other words, no matter what (t>{x) is (as long as certain boundary condi- 
tions are satisfied), the generahzed Jang equation by itself already gives a lower 
bound on the integral of the scalar curvature of g, weighted by (j). Of course the 
choice of </> affects / since </> appears in the generalized Jang equation. 

In the next couple of sections we discuss two different inequalities of the 
form 

/ Q{x)R{x)dV, (58) 

for some Q{x) > 0, where each inequality is based on one of the two proofs of 
the Riemannian Penrose inequality. The expression for Q{x) differs in the two 
cases and will be described later. However, if we then choose (j){x) = Q{x) to 
be our second equation to be coupled with the generalized Jang equation, then 
existence of such a system implies 



IGtt V IGtt 



^ > m- \/^ > I Q{x)R{x)dV ^ / Rcj) dV > Q, (59) 



proving the corresponding form of the Penrose conjecture for the original Cauchy 
data {M^,g,k). We wiU call any method of proof as above an Einstein-Hilbcrt 
action method. So far we know of only two methods of this form, the Jang- 
IMCF equations presented in the next section, and the Jang-CFM equations 
discussed in the section after that. 

Note that any inequality of the form of inequality [58] proves the Riemannian 
Penrose inequality for {M^ ,g) as a special case since then > by hypothesis. 
Thus, unless one finds a new proof of the Riemannian Penrose inequality, the 
only way to hope to prove an inequality like [58] is cither to adapt currently 
known proofs of the Riemannian Penrose inequality, as we arc about to do in 
this paper, or to use the Riemannian Penrose inequality itself. This last idea 
deserves additional consideration. 

Proof of theorem \^ Applying the divergence theorem to equation [33] gives 
us that 



/ R(l)dV > 2 / (l)q{v)dA 

= 21 <f>{h- k){v,P)dA, 
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where k is assumed to converge to zero at infinity (or have compact support), 
and df is assumed to decay at least as fast as at infinity (with reasonable 
bounds on Hess/ as well) so that 

0= lim / (t){h-k){v,v)dA, (60) 

JSr 

where v is the unit outward normal vector to E and the sphere at infinity in 



In appendix W\ we observe that 
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1/2 

{v — {v, v)v) 



and 



dA = {l + ^\dJ\^)\lY'^dA. 
Then since 1 + (jP\df\l = 1/(1 - \v\l), wc get that 



RipdV > 2 (j){h~k){v,v)dA 

M JY. 

(j){h- k){v,iy - {v,v)v) 
' (1 - bPJ^/^ 



s 



= 



if we assume that 



Sc-S7s, (1-blg)'/' 

for some smooth family of surfaces converging to E, proving the theorem. 

In the case that / is blowing up (or down) everywhere on E. then choosing 
Eg to be the level sets of / simplifies things even more since then v = ±|?;|j/. 
Then the integrand becomes ^(1 — \v\'^y^^{h — k){v, v) and equals zero with the 
usual boundary behavior since = on the boundary, \v\g is going to one on 
the boundary, and both hiy, v) and fc(i^, v) are bounded by lemma [3) Equation 
1611 is also clearly satisfied in the case of a future and past apparent horizon 
where we assume that / and (/) stay smooth and bounded, since all of the terms 
in the integrand will be bounded, and <^ = on the boundary E. 

As a final comment on theorem [S] before moving on, we fully admit that 
a better understanding of the boundary behavior of / and (j) is needed. This 
better understanding should be able to be achieved when the existence theories 
for the equations we arc proposing arc discovered. 

Before we get into applications of this theorem, it is worth noting that 



E{g,cj))^ f R^dV 
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is, up to a boundary term, the Einstein-Hilbert action of the quotiented static 
spacetime 

{S^ X M, -(jy^dt^ + g) , 

where we have turned the usual R time coordinate into an of length one to get 
a finite integral. The Einstein-Hilbert action is defined to be the total integral 
of the scalar curvature R^"^ of the spacetime. In the appendix we observe that 

and since dV^^ = cf) dV, the Einstein Hilbert action of the quotiented spacetime 
is 



/ R'^^dV^^ = / - 2A(j)) dV 

J six A/ Jm 



R(j)dV-2 / (V0,j^)gfiA. 

M JdM 

We further observe that the boundary term vanishes when (j) is harmonic on 
{M,g), as is the case in the Schwarzschild spacetime. 

Finally, the vacuum Einstein equation G = is the Euler-Lagrange equation 
which results from requiring a spacetime to be a critical point of the Einstein- 
Hilbert action. Since the Minkowski and Schwarzschild spacetimcs are the only 
vacuum static spacetimcs (with no boundary or black hole boundary) [5], it 
follows that they are the only two static spacetimcs which are critical points of 
the Einstein-Hilbert action, or equivalently E{g,(f>), since boundary terms are 
irrelevant for variations away from the boundary. 



9 The Jang-IMCF Equations 

In this section we show how inverse mean curvature flow in (M^, g) can be used 
to determine a warping factor (j) for the generalized Jang equation to get a system 
of equations which, when there are solutions, implies the Penrose conjecture for 
a single black hole when H2{M^) = 0. Alternatively, the method presented 
here has the potential to address the outermost case of the Penrose conjecture 
as stated in conjecture [5l with the additional assumption that the outermost 
generalized apparent horizon is connected. We will call the system of equations 
we are proposing in this section the Jang-IMCF equations. An important open 
problem is to flnd an existence theory for these equations. 

Before we state the Jang-IMCF equations, we need to review inverse mean 
curvature flow. As introduced by Geroch [9] and Jang-Wald [16], a smooth 
family of surfaces S](t) in {M^,g) is said to satisfy inverse mean curvature flow if 
the speed in the outward normal direction of the family of surfaces as t increases 
at each point is equal to 1/H, where H > is the mean curvature of the surface 
at that point. This flow has the important and surprising property that the 
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Hawking mass of E(f) is nondecreasing when {A{^,g) has nonnegative scalar 
curvature R. 

To be more precise, define the Hawking mass of a surface E in {M^,g) to be 



m^(S) = wgl(l--^ / H^dA 



IGtt V IStt 



where all quantities are computed in (M^,g). Then we can compute the rate 
of change of the Hawking mass of a surface when flowed out orthogonally with 
speed 77 = in {M^,g) by using the first variation formula 



dt 

the second variation formula 

d_ 

dt 

and the Gauss equation 



^{dA) = {r]H)dA = dA, 



H = -Ar] - \\UUr] - Ric{D, D)i] 



Ric(j>,P) = ii?-/? + ii/2_ l||n||2^ 

where H is the second fundamental form of S in (M'^,^), Ric is the Ricci cur- 
vature of {M^^g), and K is the Gauss curvature of E, to get 



1 1 r 2|vi?ii,^_^^, 1^, 



2 167r./s(t) H"^ ^ 2 



dt^ " " "' V 167r 
The essential assumption that E(i) is connected is used to conclude that 



/ KdA = 27rx(E(t)) < 47r 



/E(t) 

by the Gauss-Bonnet formula, which, along with 

||I|li>itr(If = i, 

allows us to conclude that 



^(mH(S(t)))>\/ME/ Adl 
dt^ JiV ^ - V 167r is(t) 167r 



S(t) 

If (A/'^, g) has nonnegative scalar curvature, then the above equation implies 
that the Hawking mass of the smooth family of surfaces determined by inverse 
mean curvature flow is nondecreasing. However, we get a more general result if 
we integrate the above equation in t and use the co-area formula 

jAdt = -dV = HdV 
V 
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to conclude that for a smooth family of surfaces satisfying inverse mean curva- 
ture flow which foliates , that 

m^(SM)-mH(S(0)) = j^H^&^ ( 1^) ^ 

where at each point x € M^, H is the mean curvature of the surface S](t) through 
the point x. Then assuming that S(0) has H = and is area outerminimizing 
and that {M^,g) is sufficiently asymptotically flat, we conclude our main result 
that 

where A is the area of S = 5](0), the area outerminimizing minimal boimdary 
of and 



Qix) = H 



\m\s 



Wit 

More generally, Huisken-Ilmanen [13] observed that there exists a weak no- 
tion of inverse mean curvature flow in which the surfaces S(t) jump outward 
to their outermost minimal area enclosures whenever they are not already that 
surface. A key step in their approach is to represent the family of surfaces 
T,{t) as the levels sets of a real- valued function u{x) on called the level set 
function. Then if 

T,(t) = d{x I u{x) < t}, 
it follows that rj = l/\Wu\g and 



H = div 



so that inverse mean curvature flow on the level sets of u{x) is equivalent to 

=\Vu\g. (63) 



\Vu\g 

Huisken-Ilmancn then proceed to define a notion of weak solutions to the above 
level set equation using an energy minimization technique. These solutions have 
"jump regions" where Vw = corresponding to where the family of surfaces T,{t) 
is not continuously varying but instead "jumps" over these regions. Further- 
more, Huisken-Ilmanen, using elliptic regularization, proved that weak solutions 
of their inverse mean curvature flow always exist. We refer the reader to their 
beautiful work [14| . However, using their generalized inverse mean curvature 
flow, wc achieve the following theorem. 

Theorem 7 Given an asymptotically flat {M^ ,g) with H2{M^) — and a min- 
imal connected boundary S which bounds an interior region, then 



^ > [ Q{x)(^)dVix) (64) 
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where A is the area of the outermost minimal area enclosure E = dU^ of E, R 
is the scalar curvature and dV is the volume form of g, and 



where u{x) is a weak solution to Huisken-Ilmanen inverse mean curvature flow 
equalling zero on E. 

Proof: As described in Huisken-Ilmanen's paper, if E is not already its own 
outermost minimal area enclosure, it immediately jumps to it. During this initial 
jump, but only on this first jump, the area of the surface may decrease. Hence, 
A must be defined to be the area of the outermost minimal area enclosure E of 
E, which also has zero mean curvature by the maximum principle using E as 
a barrier. Also, since H2{M^) = 0, it follows that S is connected. Since each 
component of E bounds a region, it follows that if E did have more than one 
connected component, all of the components except for one could be removed 
(by either filling in holes or removing disconnected regions), thereby decreasing 
the area. Then starting the flow at E, our previous calculations generalize. 
The condition that H2{M^) = is also used to guarantee that each E(f) is 
connected after each jump and therefore has Eulcr characteristic < 2 as needed 
in the computation of the rate of change of the Hawking masses of E(t). 

By the first variation formula mentioned earlier in this section, inverse mean 
curvature flow grows the area form exponentially. Hence, 



which equals the desired result by equation 1631 

Theorem [7] deserves careful consideration. In the case that R > 0, we re- 
cover a Riemannian Penrose inequality for a single black hole. More generally, 
however, since ^/167r is energy density, we see that we have a kind of integral 
of energy density on the right hand side of equation 1641 modified by the factor 
Q(x). On the flat metric on and starting inverse mean curvature flow at a 
point, Q{x) = 1, and on the Schwarzschild metric, Q{x) is the harmonic function 
going to one at infinity and equally zero on the minimal neck. This last fact, 
which can be verified by direct calculation, will turn out to be important since 
this harmonic function also equals the warping factor 0(x) in the Schwarzschild 
metric. 




(65) 



|E(t)|,=ie*. 



Thus, we have that 
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Theorem [7] and theorem IH] together motivate the system of equations, 



= trg(/i-fc) 



(66) 




(67) 



(68) 
(69) 



= Q, 




where for our later convenience we choose c = y . The first equation is the 
generahzed Jang equation again. The second equation is the level set formu- 
lation of inverse mean curvature flow on (M^,g). The third equation is the 
definition of Q{x) in terms of the inverse mean curvature flow level set function 
u{x). The new equation, then, is the fourth equation, which sets 0(a;) equal to 
Q{x), up to a constant. Then by theorems [7| and [HI we conclude that 



167r Jm 

Then recalling that g measures areas at least as large as g does, we have 
that 



Eg is the outermost minimal area enclosure of E in {M^,g). Recall also that 
since H2{M^) = 0, E connected (and bounding a region) imphes that both Eg 
and Eg are also connected (and bound a region). Hence, if we can solve the 
above system with boundary conditions so that E has zero mean curvature in 
{M^,g) and so that the total masses of {M^,g) and {M^,g) are the same, then 
we would be able to conclude that 



which would prove the Penrose conjecture for a single black hole in the case 
that H2iAP) = 0. 

In the case of equality in the above inequalities, {M^,g) has to be a time 
symmetric slice of the Schwarzschild spacetime by the original Huisken-Ilmanen 
result. Thus, inverse mean curvature flow yields precisely the spherically sym- 
metric spheres of Schwarzschild, so u is easy to compute. Direct computation 
then reveals that Q{x) is the harmonic function in {M^,g) which equals zero 
on E(0) and goes to one at infinity. Since (j) equals Q (up to a multiplica- 
tive constant, which is irrelevant after a constant rescaling of the time coor- 
dinate in what follows), wc get that (R x , —cfPdt'^ + cfj is isometric to a 
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Schwarzschild spacetime. Hence, g = g — (fy^df"^ is the induced metric on a 
slice of Schwarzschild with graph function f{x). Finally, examining the case of 
equality of theorem [5] (and that theorem's use of the generalized Schoen-Yau 
identity) forces \\h — = 0, which of course implies that fcy = hij. Hence, 
the original Cauchy data {M^,g,k) is the induced Cauchy data on a slice of a 
Schwarzschild spacetime with graph function f{x). 

Thus, understanding this system of equations, and whatever existence theory 
might be associated with it, is a very interesting and important open problem. 
A first step is to observe that Q does not need to be defined in the system. 
Hence, our system is equivalent to 

The Jang - Inverse Mean Curvature Flow Equations 



= tTg{h-k) (70) 



IVtzl, = divfj^) (71) 
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\Vu\, e"/^ (72) 



where we recall that 
and 



h = 



0Hess/ + {df ®d(j) + d(i)® df) 



(l+(/)2|rfJ|2)l/2 

which can be thought of as three equations and three free functions /, m, and (f> 
on the original Cauchy data {M^ , g, k). 

In fact, the third equation, eauationlT^ can be used to solve for in terms of 
u, du, and df . The purpose of this is to recognize that the Jang-IMCF equations 
may also be thought of as two equations and two free functions / and u once we 
substitute for (f>. Since only first derivatives of / and u appear in the expression 
for (j) below, the resulting equivalent system is two second order equations in / 
and 4>. 

Unfortunately, the expression for in terms of / and u on {M^,g) is a bit 
messy, but at least it is explicit. From equation I72i we get 

= \du\g e"/2 (73) 

which is simple enough except that (j) also appears in the expression for g. Next 
we note that 



2 _ T^lj, 



\du\g = g'^u.Uj 



l + <j>^\df\l 
^Hdf,du)l 



2 
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which, when combined with equation 1731 gives us 

It follows that 0^ solves the quadratic equation, 

\df\l-<l)^ + B-<j>^-e^\du\l^O, 
where B = 1 + {{df, du)^ - \df\l\du\l). Thus, 



-B+ JB'- +'le^\df\l\du\l 



(74) 



which is clearly always nonnegative (and where we disregard the negative square 
root in the quadratic formula since that solution is nonpositive) . Thus, an 
equivalent formulation of the Jang-IMCF equations is 



= trg{h-k) (75) 
'Vu\ 



|Vub = div( J^l (76) 



IS 

where 

.g = . 9 + 02^/2, (77) 
<j>Resfif + [df (g) d(p + d(p ®df) 



(79) 



\dfU 
and 

B = l + e^{{df,du)l-\df\l\du\l), (80) 

which can be thought of as two equations and two free functions / and u on the 
original Cauchy data {M^,g, k). (In equation [79l = \du\ge'^/^ when \df\g = 
by equation [73]) . 

We end this section with a general discussion of the some of the challenges 
involved in finding an existence theory for the Jang-IMCF equations. First, 
note that these equations reduce to the Huisken-Ilmanen IMCF equation on 
{M^,g) when trg(A;) = since then we can choose / = (which implies g = g) 
to satisfy the generalized Jang equation (equation [75]). Thus, clearly a notion of 
a weak solution to this system of equations is required. Furthermore, the notion 
of "jumps" must still be involved when there are regions in which du = 0. Note 
that when du — 0, then = 0. Thus, if / stays smooth and bounded, it would 
follow that h ~ 0, which means that the generalized Jang equation cannot be 
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solved unless ti g{k) = in this region as well. If tr g{k) 7^ in this region, then 
this would suggest that / needs to be unbounded or undefined in this region. 
Clearly this is an important issue to understand. 

Given these and other considerations, one might be tempted to be pessimistic 
about finding a general existence theory for the Jang-IMCF equations. In fact, 
it was once the case that most were pessimistic about the original inverse mean 
curvature flow proposed by Geroch [9] , right up until Huiskcn-Ilmanen |14| found 
an amazingly beautiful and natural existence theory for a generalized version of 
inverse mean curvature flow. Thus, there is also precedent for optimism. 



In this section we comment that there is at least one other Einstcin-Hilbert 
action method in addition to the Jang-IMCF equations. So far we have seen 
how the Penrose conjecture would follow from a general existence theory for 
the Jang-Zero Divergence equations presented in section [7] or, for a single black 
hole in dimension three, from a general existence theory for the Jang-IMCF 
equations presented in section[9l In this section, we briefly discuss a third system 
of equations whose existence theory would also imply the Penrose conjecture. 
The precise statement of this third system is a bit laborious and so we do not 
state it here, but only describe it and the additional considerations it involves. 
In section |51 we explained how any inequality of the form 



for some Q{x) > 0, leads to a system of equations which implies the Penrose 
conjecture. The first equation in the system is the generalized Jang equation 
and the second equation in the system is simply (j){x) = Q{x) (times a constant 
if one likes). In section [5] we pursued this approach in detail for the Huiskcn- 
Ilmanen inverse mean curvature flow. 

Bray's proof [2] of the Riemannian Penrose inequality, when revisited, also 
yields an inequality of the form of equation 1811 This proof of the Rieman- 
nian Penrose inequality involves a conformal flow of metrics (CFM) which flows 
an initial asymptotically flat metric with nonnegative scalar curvature to a 
Schwarzschild metric in the limit as the flow parameter goes to infinity. Further- 
more, the area of the horizon stays constant, and (by the Riemannian positive 
mass theorem it turns out that) the total mass is nonincreasing during the flow. 

To generalize the conformal flow of metrics (CFM) proof to get an inequality 
as in equation 1811 we first need to generalize the positive mass theorem to get 
an inequality of the form 



Jm 

for some Q{x) > on some (M, g). Witten's spinor proof [33] of the Riemannian 
positive mass theorem provides such a result, for example, at least whenever a 



10 The Jang-CFM Equations 




(81) 




(82) 
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spinor solution to the Dirac equation exists (since we are not assuming R > 
anymore, there is an issue now). Also, a result of this type can be found by 
multiplying the metric g by a conformal factor to achieve zero scalar curvature 
globally (when such a factor exists), and then measuring how much the mass 
changes. This last idea is made precise by Jauregui in [T7| . Finally, one can also 
use inverse mean curvature flow starting from any point to prove an inequality of 
the above form in cquation l82l This third approach currently has the advantage 
over the first two in that it is known to work in all cases in dimension three by 
the previous section and the work of Huisken and Ilmanen [T3] and Streets [3T] . 
In the conformal flow of metrics {M,gt) with total masses m{t), 



m'{t) - -im(i). (83) 



Hence, 



— > / -mm (84) 



since the areas A{t) of the horizons of (M, gt) stay constant and the flow of met- 
rics converges to Schwarzschild where m — -y/^^ = 0. Then plugging equation 
1821 into the above equation and accounting how the scalar curvature transforms 
conformally gives a result of the desired form in equation 1811 Hence, modulo 
the existence questions needed to get an equality of the form of equation l82i we 
get a generalization of the Ricmannian Penrose inequality. 

One difference between the Jang-CFM equations and the Jang-IMCF equa- 
tions, however, is that the Jang-CFM equations are not local. That is, Q{x) in 
this case does not satisfy a local p.d.e. at each point and instead has a more 
complicated expression. Hence, for the Jang-CFM equations to have an exis- 
tence theory, the theory would have to work for a wide range of possible Q. 
On the other hand, the Q{x) from the Jang-CFM equations has the potential 
to have better regularity than the Q{x) from the Jang-IMCF equations. For 
the Jang-IMCF equations, ccj) ~ Q is not necessarily continuous or even posi- 
tive, and in fact equals zero in jump regions of the inverse mean curvature flow 
on {M^,g), which as discussed at the end of the previous section, introduces 
additional analytical challenges. 



11 Open Problems 

The two most interesting and important open problems discussed in this paper 
are finding an existence theory for the Jang-Zcro Divergence Equations (which 
would prove the Penrose conjecture) and finding an existence theory for the 
Jang-IMCF equations (which would prove the Penrose conjecture for a single 
black hole when H2{M^) = 0). Another interesting problem is to find a general 
existence theory for any Einstein-Hilbert action method as long as the associated 
Q{x) has certain properties. If the Q{x) from the Jang-CFM equations qualified 
for such a theory, this would also prove the Penrose conjecture. 
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Another interesting problem is to find additional Einstein-Hilbert action 
methods by finding new inequalities of the form of equation 1811 Since the 
special case of _R = implies the Riemannian Penrose inequality, one would 
either have to find a new proof of the Riemannian Penrose inequality or use the 
Riemannian Penrose inequality itself to find such a generalization. There may 
be reasonable ideas to try in this latter approach. 

There is also the question of the physical interpretation of inequalities of the 
form of equation 1811 Rewriting the inequality gives us 



which could be interpreted as saying that the total mass of a time-symmetric 
slice of a spacetime (not necessarily with nonnegative energy density) is at least 
equal to the mass contributed by the black holes (the first term) plus a weighted 
integral of the energy density (the second term), since energy density at each 
point can be interpreted as /2 = R/ldn. The purpose of Q can be interpreted as 
the need to account for potential energy. Also, Q should go to zero (and does 
in the IMCF and CFM cases) at and inside the horizons of the black holes since 
matter inside black holes should not affect the total mass. 

We also believe that the generalized Schoen-Yau identity and the general- 
ized Jang equation have much potential for many possible applications in the 
study of general relativity. One point of view, for example, is that the Jang- 
Zero Divergence equations give a canonical way of embedding Cauchy data 
(M^, g, k) into a static spacetime. If one is interested in understanding how the 
initial Cauchy data evolves under the vacuum Einstein equations, or some other 
equation coupled with the Einstein equation, then one could compute how the 
canonical static metrics associated with the evolving Cauchy data slices evolve. 
One nice property of this approach is that if the initial Cauchy data is a slice 
of the Schwarzschild spacetime, and we are solving the vacuum Einstein equa- 
tions for example, then while the Cauchy data is evolving in what might appear 
to be complicated ways, the associated canonical static spacetime remains the 
Schwarzschild spacetime. Also, since the generalized Jang equation blows up 
on horizons, this method could only be used to study the exterior region of 
spacetimes outside the apparent horizons of black holes. There may be some 
advantages to this restriction if this becomes a natural way to avoid spacetime 
singularities. 

A Introduction to the Appendices 

The target audience of these appendices are graduate students and other re- 
searchers who are interested in entering geometric relativity as a field to study. 
As such, we have included more detail in these calculations than is typical. We 
justify this choice in part with the fact that there are so many computations, 
many people would have a hard time duplicating all of these computations in a 
reasonable amount of time, even with well chosen hints. We also hope that these 
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appendices will be useful to students and researchers who are interested in prac- 
ticing their computational skills. We recommend the book "Semi-Riemannian 
Geometry with Applications to Relativity" by Barrett O'Neill as an excellent 
introduction to the differential geometry of general relativity, and we mostly 
follow that book's notation here. Readers should go through the calculations in 
these appendices in order since notational conventions which are established in 
one appendix apply to the appendices which follow as well. 

The authors would like to thank Alan Parry for helping with the TeXing 
and Jeff Jauregui for helpful comments improving the readability of these ap- 
pendices. 

B Curvature of Static Spacetimes 

In this section we compute the Einstein curvature, Ricci curvature, and scalar 
curvature of the general static spacetimc metric 

on R X M^, where i £ R, a; 6 and g is a positive definite metric on M^. 

First, choose a coordinate chart on {M^, g) with coordinates {x^, a;^, x^) and 
let = t. Then {da = gf^la^o ^ basis of the tangent plane at each point 
of the spacetimc. Let {u}"}a=o be the corresponding dual basis of one forms at 
each point of the spacetimc so that ^"{dp) = (Wc use bars over these bases 
instead of tildes to be consistent with section [5] and subsequent appendices) . 
Finally, we define the components of g (which are the same as the components 
of g for tangent vectors to M^) to be 

gai3 = {da, dp)g (= g{da,dp) by convention) 

so that ,900 = -(f'ix)'^, goi = dto = 0, and g,j = {Bi, dj)g for 1 < i, j < 3. 

(Notation: We adopt the convention that Greek indices alv^rays range 
from to 3 and Latin indices always range from 1 to 3. Also, we adopt 
Einstein's convention that any time an index is both an upper and 
lower index in an expression, summation over that index is implied.) 

(Recall also that g"^ are the components of the inverse matrix of g expressed 
as a matrix at each point of the coordinate chart, and that indices of a tensor 
may be raised or lowered by contracting with g"^ or gap, respectively [23]). 

By the Koszul formula the Levi-Civita connection V of g can be ex- 
pressed in terms of its components as 

where ^ 

^ap = -jp^ {~9ae,P + gpe,a " (jaPfi) (86) 

are called the Christoffcl symbols of V. 
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(Notation: Commas denote differentiation with respect to the coor- 
dinate chart so that gaf},e = 

Plugging in our expressions for cjap^ short calculations reveal that 

^ 00 ^ ^ Oi ^ ^ iO ~ ^ ij ' 

f 0," = f ^ ^ , and that t„,i - • #, 

where cjji = ~ and, as previously stated, = g^^4>j = g^-'4>j = g^"4'a = 

4>^ since (j> is only a function of x and does not depend on t. Note that we 
are using the convention that a raised index with a bar or tilde over it denotes 
raising the index with g or g, respectively. 

Since the Lie bracket of coordinate vector fields is zero, it follows from the 
definition of the Riemann curvature tensor that 

Hence, 

For the beginner, we note that about half of the text books define the Riemann 
curvature tensor to be the negative of what we used above. However, all texts 
eventually end up with the same definition of the Ricci curvature (defined in a 
moment) which is agreed to be a positive multiple of the metric on the standard 
sphere. 

Plugging in our formulas for the Christoffcl symbols, we thus compute that 

V <^ A, ™ ' ' 

HcsSjTj^ 
from which it follow that 
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where Hess is the Hessian and A is the Laplacian on {M^,g). The second com- 
putation foUows from the first by first lowering the raised index (introducing 
a factor of —4>^), using the antisymmetry of the Riemann curvature tensor to 
switch indices, and then taking the trace of the Hessian. Wc remind the reader 
that the Latin letters j, k, m range from 1 to 3. The beginning student should 
review the definition of the Hessian, the Laplacian, and the use of normal co- 
ordinates. In this case, we note that wc may choose normal coordinates on 
such that i' j™ = at a single point. Similarly, it is straightforward to verify 
that 

-^Ojo = = -Rfejo • 

It turns out that the above components of the Riemann curvature tensor are 
all that we need to compute the Ricci curvature of the spacetime. For example, 

RiCjfc = Rajk = Rajk + ^Ojk = ^ajk + ^Gjk = RiCjfc + R^jk ■ 

The first equality is the definition of the Ricci curvature as the trace of the Rie- 
mann curvature tensor. For the second equality recall our summation convention 
for Latin and Greek indices stated above. The third equality is a consequence 
of the Gauss equation for submanifolds since the t ~ Q slice of our spacetime 
has zero second fundamental form by symmetry. The fourth equality is simply 
the definition of the Ricci curvature Ric of {M^,g). Also, 

R-icoo = ^qoo" ^ ^ooo" + ^jm = ^jm 
by antisymmetry of the Riemann curvature tensor. Finally, 

RiCjO ~ Rajo" ~ ^QjO + ^kja — 0- 

Thus, putting it all together, we have formulas for the components of the 
Ricci curvature of the static spacetime metric g = —(l){x)'^dt^ + g onR x , 
namely 

Rlcnn = 



Tj~ -o^ Hessjfe^ 
RiCjA; = RiCjk ^ — 

Ricjo = Ricoj = 

in terms of the Ricci curvature Ric of {M^,g) and the Hessian and Laplacian of 
cj) on {M\g). 

Next we can compute the scalar curvature of the spacetime by taking the 
trace of the Ricci curvature, 

i? = g^'=RiCjfe =i?-2^. 
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Hence, since G = Ric — \Rg, the components of the Einstein curvature tensor 
of the static spacetime metric are 

Goo = ^i?<^' 

Gju = R-iCjfc 
Gjo Ricoj 

as desired. 

C The Second Fundamental Form of the Graph 

In this section we will compute the second fundamental form of a space-like slice 
of the static spacetime (R x , g) , where 

(/) is a real-valued function on M, and ^ is a Riemannian metric on Af . Given a 
real- valued function / on M , define the graph map 

F : M ^Rx M (88) 

where F{x) = {f{x),x). 

As we established in section [5l we will abuse terminology slightly and always 
refer to the image of the graph map F{M) simply as M and the t ~ slice of 
the constructed spacetime as M. This notation is convenient since then {M, g) 
and (M, g) are space-like hypersurfaces of the spacetime (R x , —cjy^dt'^ + g) 
(given appropriate bounds on the gradient of /). Let tt : M i— > M be the 
projection map 7r(/(x),x) = (0,a:) to the t = slice of the spacetime. 

We repeat some notation and definitions from section [5] for clarity. Let 
Bq = dt and {di} be coordinate vectors tangent to M. Define 

9. = a, + 00 (89) 

to be the corresponding coordinate vectors tangent to M so that 7r*(9i) = di. 
Then in this coordinate chart, the components of the metrics g and g induced 
from the spacetime are 

gij^{di,d.j) and g^j ^ {di,dj) 

where the angle brackets refer to the spacetime metric. Then it follows imme- 
diately that 

9ij =9i3 ~ 4>^fifj- (90) 

We comment here that the reader should think of g, (p, and / as the variables 
that we get to choose which determine g. The metric g comes from the initial 



Hessjfc(/. fAcj) R\_ 

^^[-~2)^ 

= 
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Cauchy data (M, g, k) and (f> and / arc functions which will satisfy a system of 
equations of our choosing. 

The inverse of {gij} turns out to be 



where 



g'^=g'^+v'v\ (91) 



(92) 



(l-02|d/|2)l/2 (l + ^2|rfJ|2)l/2- 

The above computation is most easily verified at each point in normal coordi- 
nates of g at that point, where the gradient of / is assumed to lie in the first 
coordinate direction. The second part of equation [HH can be computed in the 
same manner, but where we consider that gij = gij-\-4>^ fifj and then use normal 
coordinates as before, but this time for the metric g. (For the beginner, the use 
of normal coordinates is exemplified in more detail in a moment.) 
We also define 

V = v'B^ and V = v'd^ (93) 
so that TTf{v) = V, and observe the useful identity 

{l^cb'\df\l)-{l + cf,'\df\l)^l, (94) 

which follows directly from computing the ratio of the volume forms of g and g 
two different ways, namely with respect to g and then g. 
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As established in section [51 we use the convention that a barred index (as 
in above) denotes an index raised (or lowered) by g as opposed to g. That 
is, P = g^-' fj, where as usual fj = df/dx^ in the coordinate chart. In general, 
barred quantities will be associated with the t — slice {M,g) and unbarred 
quantities will be associated with the graph slice (M, g) . 

Our next step is to compute the unit normal vector n to the graph slice 
{M,g) defined by /. It is straightforward to verify that 



do 



<f>(l^r\df\f) 



1/2 



(95) 



has unit length in the spacetime metric, is perpendicular to the tangent vectors 
di — di + fido to the graph slice, and hence must be the correct expression. 

Following our convention for the definition of the second fundamental form 
defined in equation [31 we thus have that the components of the second funda- 
mental form h are 

hij = h{di,dj) 



(Bi+fido) 



(Bi+fido) 



do 



,1/2 



(l - 
^0 + 02/^4] , d, + J, do 



dj + fjdo 



1/2 



where V is the Levi-Civita connection on our spacetime (R x M'^, g) . The third 
and fifth equalities follow from the fact that {n, dj) — Q on M. 

From the form of the above expression, we sec that the Christoffel symbols 
of the spacetime, defined and computed in appendix [b1 are going to come into 
play. From those computations, it follows that 



'-Bo 



Vs,^o = ^^o'^o 



-do 



where we remind the reader that Latin indices, when summation is implied by 
one raised and one lowered, only sum from 1 to 3, and a bar over a raised index 
indicates that the index was raised with g as opposed to g. 

It now becomes convenient to use normal coordinates on {M,g). Note that 
these are not normal coordinates on the whole spacetime, just on the t = 
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slice of the spacetime {M,g). Since by symmetry this slice has zero second 
fundamental form, Vg.dk = = at a single point of our choosing. In 

normal coordinates, derivatives of the metric components gij and g^^ are zero 
at the chosen point, thereby making the Christoffel symbols zero at that point 
as well. In addition, ~ (/m5™'^)i = fimd"^'^ = fi'' at this single arbitrary 

point. Note that fim ~ dx'dx"^ simply a coordinate chart second derivative 
in our notation. 
Hence, 

= [j + '^/^/''/'fe) do + (20(/../^ + + c^M^) dk. 

Since {do, do) — —(j)^ and {dk,dj) — gkj (which then lowers indices on other 
terms), we have computed that 

--^V, f % + + f200J^ + + <j>M^) gkj 



1/2 



0(1 "'^^M/l^) 



1/2 



1/2 
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at the chosen point in normal coordinates for (M, g). Since HesSy/ = fij as 
well at the chosen point in these normal coordinates, we have that 



, _ 0Hess,;j/ + + - (j)^df, d(t))gfifj 

at the chosen point. However, the above equation represents the components of 
the tensorial equation 



_ 0Hcss/ + d/ # + # ® df - (d/, d<^)gd/ ® df 

^ ; (98) 

by which we mean both the left and right sides of the equation arc tensors. Since 
tensorial equations may be verified in any coordinate chart, we conclude that 
equation 1981 is true at our chosen point. Since our chosen point was arbitrary. 
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it follows that equation [Ml is true at every point. Thus, equation [571 is true at 
every point as well, in any coordinate chart. 

The beginner differential geometer should take note on how using normal 
coordinates simplified these computations substantially. However, the compu- 
tation can also be done straightforwardly without using normal coordinates, 
just not as elegantly. 

As we mentioned originally in section [5l we extend h trivially in our con- 
structed static spacetime so that h{dt, ■) = 0. Note that this gives h{di,dj) = 
h{di,dj), so we can call this term hij without ambiguity. Our next goal is to 
convert the above formula for hij expressed with respect to g to one expressed 
with respect to g. 

To convert the tensor Hess/, we must recall that it is defined to be the 
covariant derivative V of the 1-tensor df . Note that df does not involve any 
metric, since df{W) = W{f) by definition. However, V does involve the metric 
g when applied to tensors. For example, in coordinates, 

H^u7=v(d/)(a„a,) 

=9,(d/(9,))-d/(Vg,a,) 

— fij ^ ^ ij ' fk 

which involves the metric g and its first derivatives by equation 1861 Hence, 

Ifeiiy / - Hess,j / = - (f,/ - r,/)/fc. (99) 

Thus, we need to compute the difference of the Christoffel symbols of g and g. 
By equations [90] and [9T1 

(aim,] +9jin,i ~gi3,m " [4''^ fifrn) ~ {<t>'^fjfm),i + {(j)^ f i f j) ,m) 
= (5''™ + u'^U™) {gim,j + 9jni.i - 9ij,m " 2(j)(j)jfifm ~ " (f^ fifmj 

-24)(j)ifjfm - (fy^fmfji - (t'^fifmi + 20(/)„i/j/j -|- (/j'^fjftm + (f^fifjni] 
= (5'=" + v'^v"') + gj^,i - g.j^ra - 2(l)^fmfzj 

-'24>4>jfifm - '^4>4>if]fm + '^4>4>jyJif]} (100) 
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so that in normal coordinates on {M,g), 



k 



( -k 'l^^\df\lf'\ 



1 - 0^ M/lg 



= /i,,i;*^ /.(^^ (101) 

(Note that the above equation is a proof of identity [3] from section [5l) Plugging 
this equation into equation [99] gives us 



Bessijf - (j)fifj{d<l), df)g = HesSy / - hijv{f). (102) 

Hence, using the above formula and equation 1941 we can transform equation 
EZlto 

, _ (^HesS;;j/ + + (j),fj - <j)'^{df, d(j))gf,fj 

~ ( ,N 1/2 

1 /2 

= (l + 02 |d/|2) (0Hess,y / - 0%i;(/) + + 0,/^) 
SO that by the definition of v in equation 1921 



1 /2 

+ 0' M/l') % =(i + M/l') (0Hcssy / + + 'i>^^j) ■ 



Thus, we sec that the second fundamental form of the graph slice expressed 
with respect to the metric g is 

- (l + 02|d/|^)V2 ' ^'"'^ 

as claimed in section [S] Technically, we've shown that the above equation is 
true at an arbitrary point in normal coordinates. Again, the above equation 
represents the components of the tensorial equation 

(/)Hess/ + {df ®d(j) + d(j)® df) 



(l + 02|d/|2)l/2 



(104) 



which is therefore true at the arbitrary point, and thus is true everywhere. 
Hence, equation 11031 is true at every point as well, in any coordinate chart. 
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D Derivation of Identities 



In this appendix wc finish the proof of the generalized Schoen-Yau identity 
sketched out in section [5] We begin witli equation |40] derived in that section, 
which we repeat for clarity: 

R = 167r(Ai- J(^;)) + (tr,/i)2-(tr,fc)2_||/^||2 + ||fc||2 
+2vitTgh) ~2v{tTgk) - 2div(/i)(i!) + 2 div(A:)(t;). 

We will convert this formula for R to an expression in terms of the g metric. 
To perform the conversion, we need several identities (originally listed in section 
[5]) for arbitrary symmetric 2-tensors k which we now prove. We continue with 
the notation established in the previous appendix, which might be thought of 
as an introduction to this appendix. 

Identity [1] 

- ||fc||^ = (tr,fc)2 - ||fc||2 + 2k{v,v)tT,k - 2\kiv, 

Proof: 

{trgkf ^ \\k\\l = {g^^h.f - g^'g^'h.kki 

= [tr, k + k{v, v)f - \\k\\l - k{v, vf - 2 \k{v, ■% 
^{tv-gkf ^\\k\\l + 2k{v,v)tv,k~2\k{v,-)\]. 

The first equality is true by definition. The second equality follows from equation 
1911 For the third equality, remember that k is defined to have zero time-time 
components and time-spatial components. Hence, k(v,w) = fc(w,w) for all w 
since v projects to v (and consequently v and v are equal except for their time 
components) . 

Identity [2] 

v(trgk) = v(trgk + k{v, v)) 

Proof: 

v{tTg k) = V (trg k + k{v, v)) 
= V {ivg k + k{v, v)) . 

The first equality was shown in the proof of identity [TJ The second equality 
follows since v and v only differ by their time components, and the function 
being differentiated does not depend on the time coordinate, by definition. 

Identity [3] 
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Proof: Sec equation IIOII 



Identity H 



div(fc)(w) = diy{k){v) + (y^k){v,v)-2\v\lk{v,^ 



+ {h{v, ■),k{v, ■))g + 2h{v, v)k{v, v) + {tvgh)k{v, v) 

Proof: 

AiY{k){v)^g'={VdMdj,v) 

— jj (h. ._r ™i- — r ""i"- "171" 

= ^{k){v) + {ti-g h)k{v, v)~(t> \df\l fc(V</., v) + {h{v, ■),k{v, ■))-, 

- (j) v{f) k{Vf,^(l)) + (^.sk){v,v) + h{v,v)k{v,v) 

- v{ff kiycj), v) + h{v, v)k{v, - v{ff k{v, V0) 
= di^(A:)(5) + (ti-g h)k{v, v) + {h{v, •), k{v, ■))§ + (^vk){v, v) 

+ 2h{v,v)k{v,v) - k (^v, {202 \df\l \v\l + 24? \df\]\ 

= (V • k){v) + (trg /i)fc(w, v) + •), ■))-g + (VbA:)(z;, v) 

+ 2/i(z;,w)fc(w,w) -2|v|^A: (^v,^ 

The first equality is the definition of divergence. The second equality is the 
definition of the covariant derivative of a tensor. The third and fourth equalities 
use Christoffel symbols as defined in the previous appendix. The sixth equality 
uses equation [5T] and then identity [31 The seventh equality is most easily seen 
by using normal coordinates with respect to g. The eighth equality combines 
terms using the fact that v is parallel to V/ in {M^ ,g) by the definition of v in 
equations 1921 and 1931 The ninth equality is the simplification 

24>' \df\l \v\l + 202 \df\l = 202 \df\l (\v\l + l) = Y^^^^ = 2 |S|2 . 



Identity [5] 
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0jDj 



Proof: First, let us clarify our notation. Recall that bars refer to the metric g. 
Hence, for example, % = gikv'' = (v, di)g, where v'' is defined in eauation[92l As 
is standard, semicolons refer to covariant differentiation (whereas commas refer 
to coordinate chart derivatives). Of course in our case, we need to specify with 
respect to which metric are we performing covariant differentiation. Hence, we 
place a bar over the semicolon to denote covariant differentiation with respect 
to g. Hence, Vi^j = {Vg.v,di)g. 

All of our computations in the proof of this identity and the two that follow 
only involve the metric so it is notationally convenient (though not really 
necessary) to use normal coordinates with respect to this metric. Then at that 
point. 



> - 0^ \df\l) 



<t>Ii3 + fi<t>] 
(l-02|d/|2)l/2 



1/2 



hW\a + ru,r){4>h) 



(1 - 02 )3/2 



<P^{df,d(f))g,fi,fj - (jy^fj 

(l-02|rfJ|2)l/2 



^'j i-i2 



1/2 



h _L|-|2 fj'^j , (P'^jdf, d4>)-gf,fj - (f>ifj + <pVifajV° 



+ \V\~- 



<i> (1-02 1^^1^)1/2 

2 -i^+,^/,(,,.)^.+ (l_^2|^^|2 



-1/2 



{wjW" d(t))gfafj - fa<P] - <PaIj) + (f^ {df , d(j))gfifj - (j)ifj] 



+ V 



,2 Vt< 



hij + Vih{v, •)j 



(t>iVj v{4)) 



((1 + \y?g)<t>^hf, - 



The above calculations follow from our formula for h in equation [551 our defini- 
tion of V in equations [52] and [Ml and the substitution 



Hi 



1 - <t>^ \dfVg 



1/2 



which we use a number of times. 
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Identity [6] 

V0 



div(fc)(u) = div(fc(t;, •)) - {h, k)g ~ {h{v, •), + ^ (^^' ' ^ 

Proof: By identity O 

div(fc(?;,-)) = div(A:)(w) + {hj,Vt^j)g 

= dd^{k){v) + {k, h)g + (fc(v, •), - ( ^ 



Identity [7] 

(Vc/c)(?;, = v{k{v, v)) - 2{h{v, ■), k{v, ■))g ~ 2h{v, v)k{v, v) + 2\v\lk [v, ^ 
Proof: By identity [5l 

so that by the definition of covariant differentiation of a symmctrie 2-tensor, 
(yyk){v,v) = v{k{v,v)) - 2k{v,Vyv) 

= v{k{v, v)) - 2{k{v, ■),h{v, ■))g - 2k{v, v)h{v, v) + 2\v\lk (^v, ^ 

proving the identity. 

Identity M 

div(fc)(w) = div(k(v, ■)) + v(k(v,v)) + k ( V, 

V / 

-{h,k)g - 2{h{v,-),kiv,-))-g + itTgh)k{v,v) 

Proof: Plugging identities [B] and [7] into identity [3] and simpUfying proves the 
identity. 

Identity [9] (The Generalized Schoen-Yau Identity) 

R = 16^^i^,-Jiv)) + \\h-k\\l + 2\q\l-^d^{(j>q) 

+ {ti-ghf - {ti'gkf + 2v{trgh - trgk) + 2k{v, v){trgh - trgfc) 

where 

q = h{v, ■) — k{v, •) = h{v, ■) ~ k{v, ■) . 
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Proof: First we recall equation [3D] derived in section [51 

R = 167r(Ai- J(^;)) + (tr,/i)2-(tr,fc)2-||/^||2 + ||fc||2 
+2vitTgh) ~2v{tTgk) - 2div(/i)(w) + 2 div(A:)(w). 

Next, we plug in identities [l] [2l and [H Note that these identities are true for 
arbitrary symmetric 2-tensors k and hence are true for h as well. Thus, 

i? = 167r(^-J(u)) 

+ {trghf - WhWl + 2h{v,v)tTgh - 2\h{v, 

-{trgkr + \\k\\l-2k{v,v)trgk + 2\k{v,-)\l 
+2v{tTgh) + 2v(h(v,v)) 
-2u(trgfc) + 2v{k{v,v)) 

-2di^(/i(u, •)) - 2v{h{v, v)) - 2h (v, ^ 



+2\\h\\l + A\h{v,-)\l~2{iVgh)h{v,v) 



+2div(fc(u, •)) + 2u(fc(v, v)) + 2k ^v, ■ 
-~2{h, k)g - A{h{v, •), k(v, ■))-9 + 2{tr-gh)k{v, v). 
Simplifying and combining terms then gives us that 
R = 167r(/i - J{v)) 

+ \\h\\l-2{h,k)g+\\k\\l 

+2\h{v,-)\l-mv,-)Mv,-))a+mv.-)\l 

-2di^(/i(z;, •)) + 2di^(fc(w, •)) - 2/1 (v, ^ ) + 2k (v, ^ 



+ {tT:ghf - {tvgkf + 2v{iVgh - trgfc) + 2k{v, v){tVgh - tr^fc) 

so that 

R = \Q^{^i~J{v)) + \\h-k\\l + 2\q\l-2'^{q)~2q{^ 

+ (trg/i)^ - (trgA:)^ + 2v{tVgh - tr^A:) + 2k{v, w)(trg/i - tr^fc) 



where 



q = h{v, ■) - k{v, •) = h{v, •) - k{v, •) 



Note that these two definitions of q exist on the entire constructed static space- 
time and are equal since both h and k are extended trivially in the time direction 
of the constructed static spacetime and v and v differ only in their time com- 
ponents. By the product rule the above equation proves the identity. 
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E The Transformation of the Mean Curvature 
of the Boundary 

In this appendix we derive the transformation formula for the mean curvature of 
the apparent horizon boundary S of (M, g) as approximated by the level sets of 
(f). Hence, to be useful, we need to assume that = on E, is strictly positive 
elsewhere, and has level sets converging smoothly to S. The discussion here 
supplements the discussion in section [6. II 
More precisely, given 

g = g + cP^df, 

we compute the mean curvature H of the level sets of (/) in (Af, g) in terms of 
the mean curvature H of those level sets in (M, g) and / and (j). 

Identity 10 Transformation of Mean Curvature Identity 

The mean curvature with respect to g = g + (f>^df^ of a level set T, of (f> is 
given by 



1/2 
1/2 

^K'^')|VE/|^(i-|r|^) 

where V is the gradient with respect to g, Vs is the gradient with respect to 
g restricted to S, H is the second fundamental form of E in {M,g) (so that 
H = trs(77)j, v is the outward unit normal vector to E in {M,g), and 

0V/ 




(l+02|V/|2)l/2 



(j)Hessf + [df ®d(t) + d(j)® df) 
(l + 02|d/|2)i/2 ■ 

Proof: To derive this identity, it is convenient to let 7 — (j)df so that g — 5+707. 
We also define 7tan = tan|, 7. Then it is a short exercise to verify that 



1 i2 \ 1/2 / ^ 
l + |7lg \ / 



(105) 



is the unit normal to E in (M, g), where 7*^ is defined to be the vector dual to 
the covector 7 with respect to g (which of course is ^V/). The above formula 
for D has the property that \v\g ~ 1 and {D,T)g ~ for all vectors T tangent 
to E. 
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To compute the mean curvature of E with respect to {M^,g) at a point 
p, choose a coordinate chart in a neighborhood of p so that the first n — 1 
coordinate directions are tangent to S. Our convention is that for a sphere in 
i?" with the standard flat metric, we choose the normal vector to be the one 
pointing outwards and the mean curvature to be positive. Hence, by identity [31 



n — 1 n — 1 n 

i,j = l ij = l 0=1 



so that 



i + i7i: 



2 \ 1/2 



i + l7i: 



7, de) +7 



|7tan|g 

7(z.)7*A 



i + i7i: 



+ l7ta 



2 \ 1/2 



1 + I7L V 1 + I7L 



Substituting 



1 + I7I 



1 + I7L 

1 + l7tan|' 



(7.rf0)g 9 
I + I7I! 



2 \ 1/2 
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we get that 



H 



1 + l7tan|^ 



1/2 



E 



i+i7ir . 



i+i7i: 



1/2 



E 



\l+ |7tan|g^ 

,2 \ 1/2 



i+i7i;; 

(7,d0) 



1 + 

,1 + |7tan|g/ 



i7 - 



n(7tan,7tan) 



f7*^ 



(rf/,7tan) 



Substituting ■y ~ (j) df, we get 



1 + K 
( - 



V 1 + I7L 



7(iy)(7,d(/))^ 
I + I7I' 



1/2 



1 



H 



tr|/i 



''n(Vs/,Vs/) \ 

1 + 0'|V/|^ j 



"0|Vs/| 



1 - 



|Vs/l: 



l + 02|V/|^ 



H<f>) 



<l>^,y{f){df,d(b), 
1 + 0^|V/|' 



(106) 



So far we have not used the assumption that E is a level set of 0, so the 
above formula would be of interest if one wanted to analyze the mean curvatures 
of a family of surfaces converging to the boundary other than the level sets of </>. 
Since in our case E is a level set of 4>, it follows that V4> \\ i' and V0 = v{(j))u. 
Hence, 



|v/i: 



Thus, if we recall that 



=^(0) 
0V/ 



'i + 0'|Vs/l: 

l + 02|V/|^ 



|v/i: 



1/2 
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and we let 



T 



tans(u)) 



1 + 0^|V/|^) 



1/2 



we get that 



( 



i + </'2|ve/i'; 



1/2 



[(i7-]i(r,T)) 



(tr|;i-/i(T,T))(i.,t;),] 




2\ 1/2 



+ 



•2K0')-|Vs/|^(i-|r|^) 



(107) 



as claimed in the identity. 

The purpose of including this identity in this paper is to help those who want 
to study the existence theories of one of the systems of equations described in 
this paper, such as the Jang - zero divergence equations. The above identity 
may be useful for understanding boundary behavior. 

In particular, to reduce the Penrose conjecture to the Riemannian Penrose 
case, it is necessary for i/ < (or something very close to this). Typically, one 
would even expect H = 0. The cases of blowup, blowdown. or bounded behavior 
everywhere are discussed in section 16.11 from a different point of view. In the 
cases of blowup or blowdown the level sets of / were used instead of those of <j>. 

However, if we are going to allow mixed blowup and blowdown behavior on 
generalized apparent horizons, then we can no longer use the level sets of / 
everywhere. However, since is always assumed to go to zero on the boundary, 
it is natural to study the mean curvatures of the level sets of (j). Since an 
existence theory is necessary before we can make very many conclusions about 
boundary behavior, we restrict this final discussion to a few observations. 

At points on the boundary in the interior of a blowup region, a blowdown 
region, or a bounded behavior region, it is plausible that v = —v, v = v, or 
V = 0, respectively, in the limit as the level sets of (j) approach the boundary. In 
those three cases, it follows that T = 0. If we further assume that the second 
term in the formula for H can be shown to be zero, we then get 



The term in brackets is then [H + (tr|, h)], [H — {ti^ h)], and [H] in those 
three respective cases. Modulo possible issues with taking limits, one could then 
use the generalized Jang equation to conclude that 



in the case of either blowup or blowdown, since g''^ = g^^ —v'^v^ . Hence, the term 
in brackets equals zero in the three respective cases of a local future apparent 
horizon, a local past apparent horizon, or a local future and past apparent 




= trg(/i ~k)^ trg(/i- fc) - (ft, - k){v,v) = tr|,(/i - k) 



(108) 
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horizon, as desired. One would then need to show that the term in front of the 
brackets still allows one to conclude H = in the limit, even when it diverges 
as the level sets of (j) approach the boundary. 

Of course the really tricky part is understanding points on the boundary 
where every open set around the point contains two or more of blowup, blow- 
down, and bounded behavior. We offer the above formula for the mean curvature 
H of the level sets of </> in case it is helpful to others who approach this problem. 
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